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Introduction 


It  is  the  purpose  of  this  report  to  sunmarize  the  electromagnetic 
research  we  have  performed  through  the  kind  support  of  the  AFOSR.  The  problems 
we  have  worked  on  pertain  to  EMP  shielding  theory  and  deal  with 

a)  the  transmission  of  electromagnetic  waves  through  small 
apertures, 

b)  the  propagation  of  electromagnetic  waves  through  chiral 
media, 

c)  the  steady  state  and  transient  electromagnetic  coupling 
through  slabs,  and 

d)  the  pondermotive  forces  and  torques  produced  by  electromagnetic 
waves. 

In  what  follows  we  describe  the  problems  we  investigated  and  the  results 
we  obtained.  The  details  of  our  work  were  circulated  among  the  cognizant  officers 
of  the  Weapons  Laboratory,  Kirtland  Air  Force  Base,  New  Mexico  and  then  published 
in  the  open  literature.  A  full  discussion  of  our  calculations  can  be  found  in 
our  published  papers,  photocopies  of  which  are  attached  to  this  report. 
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The  question  of  how  much  of  the  electromagnetic  energy  that  exists  on 
one  side  of  a  wall  can  leak  to  the  other  side  through  a  small  opening  in  the 
wall  is,  by  virtue  of  its  pratical  importance,  a  canonical  problem  in  EMP 
shielding  theory. 

For  an  opening  that  is  circular,  the  answer  is  known  from  the  celebrated 
works  of  Rayleigh,  Bethe,  Bouwkamp,  Meixner,  and  Andrejewski.  However,  for 
a  non-circular  opening  it  appears  impossible  to  solve  the  problem  exactly, 
unless  the  opening  has  a  shape  that  is  simple  enough  to  permit  a  separation 
of  the  variables  and  a  scalarization  of  the  field.  This  means  that  for  most 
openings  of  practical  interest  an  exact  solution  can  not  be  found  and  one  must 
handle  the  problem  by  a  method  of  approximation. 

One  can  formulate  the  problem  so  that  upper  and  lower  bounds  on  the  true 
(exact)  solution  and  not  the  true  solution  itself  need  be  found.  Such  a  formu¬ 
lation  can  be  based  on  Levine  and  Schwinger's  discovery  that  when  the  aperture 
(of  any  shape)  is  electrically  small  there  exists  a  variational  principle  for 
the  upper  bound  and  a  variational  principle  for  the  lower  bound.  This  varia¬ 
tional  approach  can  yield  accurate  results  but  in  general  does  not  lend  itself 
to  easy  calculation. 

In  our  work  we  discovered  a  simpler  method  of  sandwiching  the  true  solution 
between  upper  and  lower  bounds,  namely,  the  method  of  synmetrization  which  has 
yielded  important  results  in  geometry.  We  developed  the  method  in  the  direction 
of  electromagnetics  and  applied  it  to  our  aperture  problem.  One  of  our  principal 
results  is  given  by  the  simple  formula 

T<  samii 
~  27tt  (A/X^ ) 


for  the  transmission  through  an  aperture  of  area  A  and  perimeter  P.  Here  t  is 
the  transmission  coefficient  and  X  is  the  wavelength. 


From  this  formula  we  see  that  when  an  aperture  of  given  area  is 
symmetrized  the  aperture's  perimeter  decreases  and  the  maximum  possible 
transmission  decreases  as  the  shape  of  the  aperture  approaches  that  of  a 
circle.  For  example,  the  maximum  possible  transmission  decreases  as  the  shape 
of  the  aperture  is  changed  from  that  of  an  equilateral  triangle  to  that  of  a 
square  and  finally  to  that  of  a  circle. 

For  a  full  disc-*ion  of  the  method  of  symmetrization  and  its  application 
to  EMP  shielding  the  reader  is  invited  to  read:  D.  L.  Jaggard  and  C.  H.  Papas: 

"On  the  Application  of  Symmetrization  to  the  Transmission  of  Electromagnetic 
Waves  Through  Small  Convex  Apertures  of  Arbitrary  Shape",  Applied  Physics  15, 
21-25  (1978). 

Chiral  Media 

Chirality  is  a  purely  geometric  notion  which  refers  to  the  lack  of 
synmetry  of  an  object.  By  definition,  an  object  is  chiral  if  it  cannot  be 
brought  into  congruence  with  its  mirror  image  by  translation  and  rotation.  An 
object  that  is  not  chiral  is  said  to  be  achiral.  Some  chiral  objects  occur 
naturally  in  two  versions  one  being  the  chiral  object  and  the  other  being  its 
mirror  image.  Such  objects  are  said  to  be  enantiomorphs  of  each  other. 

A  chiral  object  has  the  property  of  handedness;  it  must  be  either  left- 
handed  or  right-handed.  The  handedness  of  helices  was  demonstrated  by 
W.  H.  Pickering.  His  experimental  results  showed  that  a  collection  of  randomly 
oriented  left-handed  helices  would  rotate  the  plane  of  polarization  of  a  linearly 
polarized  microwave  one  way  whereas  a  collection  of  r’ght-handed  helices  would 
rotate  the  plane  of  polarization  the  opposite  way. 

As  a  generalization  of  Pickering's  results,  our  research  led  us  to  the 
following  conjecture:  Any  medium  composed  of  randomly  oriented  equivalent 
chiral  objects  will  rotate  the  plane  of  polarization  one  way,  say,  to  the  left, 
while  a  medium  composed  of  the  enantiomorphs  of  these  objects  will  rotate  the 
plane  of  polarization  the  opposite  way,  i.e.  to  the  right.  By  examining  the 


wire  helix  and  the  wire  braid  as  chiral  objects,  we  obtained  results  which 
support  this  conjecture.  Moreover,  we  showed  that  the  pondermotive  forces  on 
the  helix  and  the  braid  tend  to  reduce  their  chirality. 

The  practical  applicability  of  chirality  to  the  measurement  of  EMP  was 
studied  and  several  experiments  based  on  the  concept  of  chirality  were  designed 

For  further  details  please  read: 

D.  L.  Jaggard,  A.  R.  Mickelson,  and  C.  H.  Papas,  "On  Electromagnetic  Waves  in 
Chiral  Media",  Applied  Physics  18,  211-216  (1979). 

Slab  Coupling 

We  studied  the  problem  of  transmission  of  steady-state  and  transient 
electromagnetic  waves  through  a  slab.  An  analytical  solution  was  obtained  for 
the  case  of  a  linear  hemogeneous,  isotropic,  highly  conducting,  infinite  slab 
excited  by  col  linear  electric  or  magnetic  dipoles. 

We  found  that  the  transmitted  components  of  the  field  are  given  as  the 
product  (steady-state  case)  or  the  convolution  (transient  case)  of  the  corre¬ 
sponding  incident  field  components  and  a  two-term  factor.  In  the  frequency 
domain  the  first  term  of  this  factor  is  exactly  the  transmission  coefficient  of 
a  plane  wave  normally  incident  on  the  slab.  The  second  term  of  this  factor 
brings  into  play  the  finite  distance  between  the  transmitting  and  receiving 
antennas  and  becomes  significant  only  when  this  distance  is  of  the  order  of,  or 
smaller  than,  the  free-space  wavelength  (steady-state  case)  or  the  spatial  wave 
length  of  the  incident  pulse  (transient  case).  Accordingly,  we  showed  that  it 
is  possible  to  obtain  plane  wave  excitation  results  even  when  the  sources  and 
the  receivers  are  located  at  finite  distances. 

We  derived  the  conditions  under  which  measurements  made  with  source  and 
receiver  at  finite  distances  are  equivalent  to  the  same  measurements  made  with 


plane  wave  excitation. 

This  work  has  been  applied  to  the  laboratory  testing  of  EMP  slab  coupling. 


For  further  details  we  refer  the  reader  to: 

6.  Franceschetti  and  C.  H.  Papas,  “Steady  State  and  Transient  Electromagnetic 
Coupling  Through  Slabs",  IEEE  Trans.  vol.AP-27,  No.  5,  September  1979. 

Pondermotive  Forces 

It  is  well  known  that  an  electromagnetic  wave  exerts  forces  and  (possibly) 
torques  on  charges  and  currents  and  that  the  mathematical  description  of  these 
mechanical  actions  is  given  by  the  Lorentz  force  equation,  which  connects  elec¬ 
tromagnetic  and  mechanical  phenomena.  If  an  electromagnetic  wave  interacts  with, 
say,  a  metal  body,  charges  and  currents  are  induced,  and  hence  the  body  becomes 
subject  to  local  forces  and  torques.  When  a  spatial  summation  is  made,  the  total 
time-dependent  (body)  force  and  (body)  torque  acting  on  the  object  are  obtained. 

It  is  usually  believed  that  these  forces  and  torques  are  "small"  and  hence 
negligible  from  the  viewpoint  of  engineering  application.  However,  this  is  not 
always  so.  As  a  matter  of  fact,  a  strong  interest  in  this  area  is  now  emerging. 

Electromagentic  pulses  (EMP)  produced  by  nuclear  explosion  in  the  upper 
atmosphere  is  a  case  of  recent  interest.  Also  satellite  applications  have 
attracted  considerable  attention,  where  the  (relatively)  small  mass  of  the  object 
and  the  cumulative  nature  of  the  mechanical  actions  (translation  or  rotation)  can 
produce  long-range  macroscopic  effects.  Mechanical  effects  can  be  used,  for 
example,  to  steer  the  satellite  by  radiating  suitably  polarized  electromagnetic 
fields  from  the  satellite.  And  even  key  operations,  as  deploying  antennas,  can 
be  performed,  in  principle,  by  injecting  transient  currents  into  the  structures. 
There  are  cases  in  which  the  movement  of  a  metal  object  in  an  electromagnetic 
field  can  produce  additional  mechanical  stresses.  For  instance,  the  propellers 
of  a  helicopter,  which  may  be  roughly  modeled  as  rotating  dipoles,  can  be  further 
stressed  when  exposed  to  an  incident  electromagnetic  field.  Clearly  related  is 
the  problem  of  the  electromagnetic  gun. 

Accordingly,  we  examined  a  number  of  very  simple  cases  of  interaction  between 
metal  objects  and  electromagnetic  fields,  and  computed  the  resulting  mechanical 


forces  and  torques.  As  the  metal  object,  we  considered  the  simplest  one  --  the 

tuned  electric  dipole  --  in  several  configurations:  stationary,  translational, 

rotational,  receiving,  transmitting,  where  the  electric  dipole  was  isolated  or 

coupled  with  a  magnetic  dipole  as  a  turnstile  antenna  or  a  chiral  antenna.  We 

computed  the  body  force  and  torque  acting  on  the  dipole  and  its  time  average, 

for  the  case  of  an  incident,  or  a  transmitted,  harmonic  electromagnetic  field. 

As  a  general  result,  the  forces  are  always  proportional  to  the  Poynting  vector 

(divided  by  the  velocity  of  light)  times  an  "equivalent  area",  which  is  four 

2 

times  larger  than  the  usual  effective  area  —  3X  /4-m  --  of  the  (tuned)  dipole. 
This  clearly  implies  that  the  mechanical  action  depends  on  the  power  scattered 
by  the  dipole.  Similarly,  the  torque  is  proportional  to  the  Poynting  vector 
(divided  by  the  angular  frequency  of  the  wave)  times  the  same  equivalent  area. 
The  existence  of  the  torque,  however,  is  not  dependent  only  on  the  symmetry 
properties  of  the  scattering  object  or  on  the  polarization  of  the  incident 
wave.  We  examined  cases  in  which  the  incident  field  —  circularly  polarized  — 
carries  an  angular  momentum  and  yet  exerts  no  torque  on  the  dipole.  Conversely, 
we  examined  cases  in  which  the  dipole  --  a  simple  short  wire  —  radiates  as  a 
scattered  field  a  linearly  polarized  wave  (in  the  far-field)  and  has  torques 
exerted  on  it.  The  full  analysis  of  a  rigorous  boundary  value  problem  --  a 
metal  sphere  in  an  incident  field  —  showed  that  the  transient  force  has  not 
always  the  same  sense  as  the  incident  radiation;  in  other  words,  the  sphere  is 
alternately  pushed  and  pulled  in  the  direction  of  the  incident  Poynting  vector. 
The  large  variety  of  effects  we  discovered  in  this  "case"  analysis  seems  to 
indicate  that  this  is  an  underdeveloped  field  of  very  interesting  research,  and 
a  number  of  general  theorems  could  probably  be  discovered,  thus  bringing  to 
a  better  understanding  of  a  wrongly  neglected  field. 

For  a  full  discussion  of  our  work  on  this  subject  please  see: 


G.  Franceschetti  and  C.  H.  Papas,  "Mechanical  Forces  and  Torques  Associated 
With  Electromagnetic  Waves",  Applied  Physics  23,  153-161  (1980). 
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Abstract.  The  transmission  of  an  electromagnetic  wave  through  a  small  aperture  in  a 
perfectly  conducting  screen  is  examined  from  the  viewpoint  of  symmetrization. 

PACS:  02.  41 


The  question  of  how  much  of  the  electromagnetic 
energy  that  exists  on  one  side  of  a  wall  can  leak  to  the 
other  side  through  a  small  opening  in  the  wall  has 
become,  by  virtue  of  its  practical  importance,  a  canoni¬ 
cal  problem  in  the  theory  of  EMP  I electromagnetic 
pulse)  interactions  [1], 

As  is  well  known,  the  earliest  calculation  of  the 
transmission  of  an  electromagnetic  wave  through  a 
small  circular  aperture  in  a  plane  screen  of  perfect 
conductivity  and  zero  thickness  was  performed  by  Lord 
Rayleigh.  Using  potential  theory,  he  calculated  the 
transmitted  field  of  a  plane  harmonic  wave  normally 
incident  on  an  electrically  small  circular  aperture  [2]. 
Years  later  Bethe  derived  expressions  for  the  polariza¬ 
bilities  and  effective  dipole  moments  of  small  circular 
apertures.  His  results  give  the  transmitted  far  field  for 
any  angle  of  incidence  but  not  the  transmitted  near  field 
[3],  Most  recently  Bouwkamp  [4],  Meixner  and 
Andrejewski  [5].  and  Andrejewski  [6]  found  an  exact 
solution  for  both  the  near  and  far  transmitted  fields  of  a 
plane  wave  normally  incident  on  a  circular  aperture. 
Aperture  problems  can.  at  least  in  principle,  be  solved 
numerically,  but  they  cannot  be  solved  analytically 
unless  the  shape  of  the  aperture  happens  to  be  simple 
enough  to  permit  a  separation  of  the  variables  and  a 
scalarization  of  the  electromagnetic  field.  However, 
from  this  it  should  not  be  inferred  that  if  the  aperture 
problem  cannot  be  solved  analytically,  a  numerical 
method  is  the  only  way  to  obtain  a  solution.  Actually,  as 
a  preferable  alternative,  one  can  reformulate  the  pro¬ 
blem  so  that  upper  and  lower  bounds  on  the  true 
solution  and  not  the  true  solution  itself  would  have  to 


be  sought.  Such  a  reformulation  can  be  based  on  Levine 
and  Schwinger's  result  that  when  the  aperture  is 
electrically  small  there  is  a  variational  principle  for  the 
upper  bound  and  another  variational  principle  for  the 
lower  bound  [7.8].  However,  this  variational  ap¬ 
proach.  which  was  used  by  Fikhmanas  and  Fridberg  to 
find  bounds  on  the  electric  and  magnetic  polarizabi¬ 
lities  of  electrically  small  apertures  [9],  does  not  lend 
itself  to  very  easy  calculation.  Accordingly,  it  is  of  some 
interest  to  try  a  simpler  method  of  sandwiching  the  true 
solution  between  upper  and  lower  bounds. 

In  this  paper  we  shall  examine  how  symmetrization. 
which  has  yielded  interesting  results  in  geometry  and 
mathematical  physics  [10],  may  be  used  to  establish 
two-sided  bounds  on  the  electric  and  magnetic  polari¬ 
zabilities  of  differently  shaped  convex  apertures  and 
thereby  estimate  their  transmission  properties  in  a 
simple  economical  manner. 

1.  Symmetrization 

Of  the  several  kinds  of  symmetrization  that  have  been 
invented  we  shall  restrict  our  attention  to  the  sym¬ 
metrization  of  a  plane  figure  with  respect  to  a  straight 
line.  To  symmetrize  a  plane  figure  with  respect  to  a 
straight  line  L.  we  suppose  the  figure  to  consist  of  line 
segments  that  are  parallel  to  each  other  and  per¬ 
pendicular  to  L  (see  Fig.  1).  We  then  shift  each  line 
segment  along  its  own  line  until  the  line  segment  is 
bisected  by  L.  The  shifted  line  segments  compose  the 
symmetrized  figure.  For  example,  a  semicircle  of  radius 
R.  when  symmetrized  with  respect  to  its  bounding 
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SYMMETRIZATION 


Fig  1.  Example  of  symmetrization  of  a  plane  figure  with  respect  to  a 
line  L  The  semi-circle  of  radius  R  ts  symmetrized  with  respect  to  its 
bounding  diameter  to  produce  an  ellipse  with  semi-axes  R  and  K/2. 
The  ellipse,  when  symmetrized,  becomes  a  circle  of  radius  R!2.  The 
area  of  each  figure  remains  constant  but  the  perimeter  decreases  with 
each  symmetrization 


diameter,  changes  into  an  ellipse  with  semiaxes  R  and 
R  2.  A  further  symmetrization  can  transform  the  ellipse 
into  a  circle  of  radius  R  2.  Symmetrization  leaves  the 
figure's  area  A  unchanged  and  decreases,  or,  more 
accurately,  never  increases  its  perimeter  P.  For  the  case 
shown  in  Fig.  1.  the  area  is  always  nR2  2  and  the 
perimeter  varies  from  (2-mt)R  for  the  semicircle  to 
2'  :tt R  for  the  circle. 

As  an  instructive  example,  we  apply  the  principle  of 
symmetrization  to  the  calculation  of  capacitance  C.  It  is 
known  that  the  symmetrization  of  a  plane  conducting 
plate  decreases  (i.e..  never  increases)  the  electrostatic 
capacity  of  the  plate  [10].  A  plane  figure  symmetrized 
infinitely  many  times  becomes  a  circle  and,  con¬ 
sequently.  of  all  conducting  plates  of  a  given  area  the 
circular  plate  has  the  minimum  capacity.  Accordingly. 

CgC,„.  (II 

w  here  C  denotes  the  electrostatic  capacitance  of  a  plane 
conducting  plate  and  C,„  denotes  the  electrostatic 
capacitance  of  the  circular  plate  of  radius  rin.  that  has 
been  obtained  by  completely  symmetrizing  the  original 
plate  This  places  a  lower  bound  on  C.  To  obtain  an 
upper  bound,  we  invoke  the  conjecture  that  of  all  plates 
with  a  given  perimeter,  the  circular  plate  has  the 
maximum  capacitance  [10],  Thus  we  find 

C„„,gC.  (2) 

where  C„u,  is  the  electrostatic  capacitance  of  a  circular 
plate  of  radius  rOUI.  whose  perimeter  is  equal  to  that  of 


the  perimeter  of  the  original  plate  From  ( 1 1  and  (2i  it 
follows  that 

C„$CSCM.  (3 1 

Since  we  have 

r,„  =  (A  7t|’  *  ( 4 1 

^  ‘5> 

and  the  electrostatic  capacitance  of  a  circular  plate 

(disk)  in  MK.S  units  is  given  by 

C  =  8e0a.  (oi 

where  a  is  the  radius  of  the  disk  and  r,0  is  the  dielectric 
constant  of  free  space,  upon  replacing  a  by  rin  and  roul 
we  obtain  from  (3l-(6)  that  the  capacitance  C  of  a  plate 
of  area  A  and  perimeter  P  is  delimited  by 

(A  n)' 2  §  C/Se,,  S  P  2n.  (7| 

Here  e0  =  (367t)' 1  x  10“ 9  farads  per  meter. 

Both  Maxwell  [11]  and  Rayleigh  [12]  made  unproven 
statements  concerning  bounds  on  the  capacitance  of 
plates,  which  agree  with  (7).  Moreover,  the  capacitance 
of  an  elliptic  plate  of  eccentricity  e.  as  given  by 

^ellipse  K 

=  (.47i|‘  2(1  -e2)~'  *,2K[e2)-^r~(A  rt)1  :(1  ~  e:  b4i. 

(81 

where  K(e2)  is  the  complete  elliptic  integral  of  the  first 
kind  [12],  clearly  satisfies  the  left  side  of  ( 7 1.  To  show 
that  it  also  satisfies  the  right  side  we  only  need  to  recall 
that  for  an  ellipse 

^ellipse  ~n 

=2(,4-7t)1 2E(e2)(l-e2r14  rt-— (.4  rrF  2 ( 1  -t-3c2  64) . 

(9i 

where  £(e2)  is  the  complete  elliptic  integral  of  the 
second  kind. 

By  virtue  of  the  apparent  validity  of  ( 7 )  for  the 
capacitance  of  plates  of  arbitrary  size  and  shape  we  are 
led  to  believe  that  other  quantities  of  physical  interest 
may  be  similarly  sandwiched  between  bounds  involving 
only  the  purely  geometric  parameters  .4  and  P. 

2.  Polarizabilities  and  Transmission  Coefficients 
of  Small  Apertures 

Let  us  now  consider  the  transmission  of  electromag¬ 
netic  energy  through  an  electrically  small  aperture 
which  is  located  in  a  plane  screen  of  perfect  con¬ 
ductivity  and  zero  thickness.  Since  the  aperture  is  small, 
the  fields  on  the  shadow  side  of  the  screen  appear  to 
emanate  from  dipoles  located  in  the  aperture.  These 
electric  and  magnetic  dipoles,  having  moments  p  and  m. 
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radiate  in  free  space  and  are  linearly  related  to  the 
incident  traveling  wave  through  the  vector  electric 
polarizability  with  components  a,  and  the  dyadic 
magnetic  polarizability  with  components  f},r  That  is, 

p,  =  e„x,Er  (i=l.  2, 3)  (10) 

I  C'=1.2)  (ID 

i 

where  p0  =  47tx  10' 7  henries  per  meter.  The  incident 
electric  and  magnetic  fields  are  plane  waves  of  the  form 
exp[i(k  r-  wt )]  and  H'"c  exp[i(k  r  — tor)],  where  r 
is  the  position  vector,  k  is  the  wave  vector  and  u>  is  the 
frequency  We  emphasize  here  that  "electrically  small" 
means  that  kL  <S  I  where  L  is  the  maximum  linear 
dimension  of  the  aperture. 

For  a  circular  aperture  of  radius  a  the  polarizabilities 
are  given  by  the  simple  expressions 

0=1. 2.3)  (12) 

(!,./=  1.2)  (13) 


The  values  1.  2.  and  3  correspond  respectively  to  the 
directions  e,,.  e  .  and  e„.  The  aperture  plane  is  defined  by 
the  unit  vectors  ev  and  e  and  the  normal  (pointing 
toward  the  shadow  sidel  is  defined  by  e„  =  e,  xe  (see 
Fig  2 1  The  polarizabilities  are  defined  here  for  incident 
traveling  waves  and  for  dipoles  radiating  in  free  space. 
For  short  circuit  incident  fields  and  for  dipoles  radiat¬ 
ing  in  the  presence  of  a  conducting  wall,  all  values  of  the 
polarizabilities  should  be  divided  by  the  numeric  4. 


a',rclc  =  faJhl3 
where 


1  ,  =j 

0  i*i. 


j.  □' 


za‘Q^ 


F Z  l  nil  vectors  <  and  i  lie  in  the  aperture  plane,  and  r„.  For 
polarization  H’"‘  is  always  parallel  lo  the  aperture  plane  and  makes 
angle  /  with  respect  toi  tor  i  polarization  E'"‘  is  always  parallel  to 
the  aperture  plane  and  makes  angle  /  with  respect  to  ( 


For  elliptic  apertures  with  semi-axes  a  and  b  along  i. 
and  e±,  respectively,  we  have 


ab* 

*  ~  3  £<e2) * 


( 1 4 1 


K 


ipse _ 


4rt  ab2e2 

3  (1 -e2)[K(e2)-£(e2)]  “ 

47t  ab2e 2 

3  £(e2)-(l -e2)K(e2)  ,2‘ 


(15) 

(161 


where  e  =  (l  -b2  a2)1  2  is  the  eccentricity  of  the  ellipse 
and  A.'(e2)  and  £(e2)  are  elliptic  integrals  of  the  first  and 
second  kind  [13]. 

The  transmission  coefficient  r  is  defined  as  the  ratio  of 
the  total  far-field  power  transmitted  through  the 
aperture  divided  by  the  total  power  incident  on  the 
aperture.  For  the  case  where  the  principal  axes  of 
magnetic  polarizability  dyadic  correspond  to  e  and 
e  we  find 


k* 


12n.4 


35  sin*  y 


-K/fj,  sin2y +  ^;:cos2y)|C°j  ;jl  ( 1 7 1 

for  j1)  polarization  [14],  Here  y  is  the  angle  of 

incidence,  i.e.  the  angle  between  k  and  e„.  and  y  is  the 
angle  between  H'"‘  and  e_  for  parallel  polarization  and 
is  the  angle  between  £mc  and  eu  for  perpendicular 
polarization  (see  Fig.  2). 


3.  Bounds  on  Polarizabilities 
and  Transmission  Coefficient 


Imitating  the  procedure  we  followed  to  establish 
bounds  on  the  capacitance  of  plates,  we  now  construct 
bounds  on  the  mean  magnetic  polarizability  jim  of  a 
convex  aperture  by  replacing  the  radius  a.  which 
appears  in  expression  (13)  for  the  polarizability  of  a 
circular  aperture,  by  rm  14)  and  foul  (5)  of  the  aperture 
Thus  we  are  led  to  the  conjecture 


16,  A 


(181 


where  by  definition  /?„  =  (/?,,  +  /J2.)  2. 

To  test  the  plausibility  of  (18)  we  examine  several 
special  cases.  For  the  elliptic  aperture  of  small  eccen¬ 
tricity  (c«  I).  (18)  becomes 

(15)  and  116)  vield 


/<„  = 


16, .4 


(;l  (' 


32 


1 20 1 
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and  thus  we  clearly  see  that  ( 18)  is  satisfied  in  the  case  of 
mildly  eccentric  ellipses,  ft  can  also  be  shown  that  (18) 
holds  true  for  elliptic  apertures  of  arbitrary  eccentricity 
(Ogeg  1)  and  for  other  convex  apertures  such  as  the 
rectangular  and  the  rhombical  aperture  [14, 15],  The 
fact  that  these  test  cases  are  in  complete  agreement  with 
(18)  leads  us  to  believe  that  the  assertion  (18)  is  valid  for 
all  convex  apertures. 

Accepting  the  general  validity  of  (18)  and  recalling  that 
symmetrization  reduces  P  without  changing  A  we 
conclude  that  of  all  convex  apertures  of  fixed  area  A  the 
circular  aperture  possesses  the  smallest  mean  magnetic 
polarizability. 

The  electric  polarizability  contributes  to  transmission 
through  small  apertures  only  when  the  incident  wave  is 
obliquely  incident  and  polarized  parallel  to  the  plane  of 
incidence.  To  construct  bounds  for  the  electric  polariza¬ 
bility  we  note  that,  for  a  circular  aperture  of  radius  a 
and  area  A.  (12)  can  be  written  as 


c,,c|«=  8  A; 

3ir  a 


(21) 


Then  by  replacing  the  radius  a  of  this  expression  by 
rw  l4i  and  ruut  1 5 1  we  arrive  at  the  conjecture 


16.4-  J,( 

- <  2  ,  <  -  — 

3  it  P  3  '  tr  I 


(221 


To  test  the  plausibilitv  of(22l  we  again  consider  the  case 
of  a  mildly  eccentric  ellipse  (e«ll.  In  this  Case  (22) 
becomes 


3'  n 


<  <  _ 
=  3  =  3 


(23) 


and  from  ( 14i  we  have 


~  ellipse  _ 


8;  4 


1  - 


(24) 


Obviously,  expression  (24)  is  equal  to  the  lower  bound 
in  (23i  Furthermore,  with  the  aid  of  (14)  it  can  be 
verified  that  the  lower  bound  in  (22)  is  precisely  the 
v  alue  o  f  the  electric  polarizability  of  ellipses  of  arbitrary 
eccentricity  [9,15].  Also  we  note  that  the  electrical 
polarizabilities  of  rectangular  and  rhombical  apertures 
satisfy  (22)  [14. 15], 

Assuming  the  validity  of  (22)  and  invoking  sym¬ 
metrization.  we  find  that  of  all  convex  apertures  of  fixed 
area  the  circular  aperture  possesses  the  largest  electric 
polarizability 

The  bounds  that  have  been  proposed  for  the  electric 
( 22 1  and  mean  magnetic  ( 1 8 1  polarizabilities  can  be  used 
to  obtain  bounds  on  the  transmission  coefficient  1 17).  In 
some  modern  applications  the  quantity  of  interest  is  the 
upper  bound  for  the  case  where  the  incident  wave  is 


directed  and  polarized  to  maximize  the  transmission 
through  the  given  aperture.  Clearly,  maximum  possible 
transmission  through  a  given  aperture  occurs  when  the 
incident  wave  is  parallel  polarized  and  is  made  to  fall  on 
the  aperture  at  grazing  incidence.  To  find  the  upper 
bound  for  maximum  possible  transmission  we  use  (22) 
and  note  that  rCill  5  rin.  Thus 


z  ■  2  ^64(A\3  ^M(P\b 
xtsin  ,■<—  —  ^  -T-  r—  ■ 

3  9\nj  ~  9  \2  n) 


(25) 


Moreover,  in  view  of  (18)  we  can  write 


sin-y  +  ^;:cos2yS 


1024/  P  \6 
9  \2r  / 


(26) 


Substituting  (25)  and  (26)  into  ( 1 7 )  we  thus  obtain  the 
following  expression  for  the  maximum  possible  trans¬ 
mission  through  a  small  aperture  of  area  .4  and 
perimeter  P 


■'  68(/>  /  )b 
T=  27ff3|.4  /.-) ' 


(27) 


where  /.  =  2rr  k  is  the  wavelength  of  the  incident 
radiation 

Since  symmetrization  reduces  P  and  keeps  A  un¬ 
changed  we  see  from  (27)  that  the  maximum  possible 
transmission  decreases  as  the  aperture  is  symmetrized 
That  is.  the  maximum  possible  transmission  decreases 
as  the  shape  of  the  aperture  approaches  that  of  a  circle 
[16] 


C  .inclusions 

We  suggest  bounds  for  the  polarizabilities  of  a  small 
convex  aperture  of  arbitrary  shape  and  given  area  and 
thus  find  plausible  upper  and  lower  bounds  on  its 
transmission  coefficient.  Symmetrizing  the  aperture  we 
see  that  the  aperture's  perimeter  decreases  and  accord¬ 
ingly  the  maximum  possible  transmission  decreases  as 
the  shape  of  the  aperture  approaches  that  of  a  circle. 
For  example,  the  maximum  possible  transmission  de¬ 
creases  as  the  shape  of  the  aperture  is  changed  from  that 
of  an  equilateral  triangle  to  that  of  a  square  and  finally 
to  that  of  a  circle. 

The  bounds  are  simple  to  evaluate  from  a  know  ledge  of 
the  aperture's  area  and  perimeter  and  therein  lies  the 
desirability  and  economy  of  this  method. 

It  appears  that  this  method  of  estimation  can  be 
generalized  to  handle  other  boundary-value  problems 
and  thus  provide  information  as  to  how  their  solutions 
are  modified  when  there  is  a  change  of  shape 
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Abstract.  We  analyze  the  propagation  of  electromagnetic  waves  through  chiral  media,  i.e.. 
through  composite  media  consisting  of  macroscopic  chiral  objects  random!)  embedded  in  a 
dielectric  The  peculiar  effects  that  such  media  have  on  the  polarization  properties  of  the 
waves  are  placed  in  evidence  To  demonstrate  ihe  physical  basis  of  these  effects,  a  specific 
example,  chosen  for  its  analytical  simplicity,  is  worked  out  from  first  principles 

PACS:  02.  41 


Chirality  is  quite  common.  It  occurs  not  only  in  nature 
but  also  in  works  of  an  and  architecture  as  well  as  in 
manufactured  articles  [1.2].  In  nature  we  find  chirality 
on  a  molar  scale  in.  for  example,  snails,  flowers,  and 
vines,  and  on  a  molecular  scale  in  such  substances  as 
grape  sugar  and  fruit  sugar  Moreover,  chiraht)  is  an 
operational  feature  of  such  manufactured  articles  as 
screws,  springs,  and  golf  clubs. 

Since  chiralit)  begets  handedness  and  handedness 
begets  optical  activity,  it  is  not  surprising  that  the 
interaction  between  an  electromagnetic  wave  and  a 
collection  of  randomly  oriented  chiral  objects  can  be 
such  as  to  rotate  the  plane  of  polarization  of  the  wave 
to  the  right  or  to  the  left  depending  on  the  handedness 
of  the  objects. 

The  concept  of  chirality  is  not  new.  nor  has  it  been  ig¬ 
nored  Since  the  early  part  of  the  nineteenth  century,  it 
has  played  an  increasingly  important  role  in  chemistry 
[3-5],  optics  [6.7],  and  elementary  particle  physics 
[8].  In  1811  Arago  [9]  discovered  that  crystals  of 
quartz  rotate  the  plane  of  polarization  of  plane  polar¬ 
ized  light  and  hence  are  optically  active.  Short!) 
thereafter,  circa  1815.  Biot  [10]  discovered  that  this 
optical  activity  is  not  restricted  to  crystalline  solids  but 
appears  as  well  in  other  media  such  as  oil  of  turpentine 
and  aqueous  solutions  of  tartaric  acid.  These  discov¬ 
eries  led  to  the  fundamental  problem  of  determining 
the  basic  cause  of  optical  activity  In  1848  Louis 
Pasteur  [3]  solved  the  problem  by  postulating  that  the 
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optical  activity  of  a  medium  is  caused  by  the  chirality 
of  its  molecules  Thus.  Pasteur  introduced  geometry 
into  chemistry  and  originated  the  branch  of  chemistry 
we  now  call  stereochemistry  More  recently,  m  1920 
and  1922.  Lindman  [11.  12]  devised  a  macroscopic 
Imolarl  model  for  the  phenomenon  by  using  mi¬ 
crowaves  instead  t.  light,  and  wire  spirals  instead  of 
chiral  molecules  The  validity  of  the  model  was  verified 
a  few  years  later  by  Pickering  [13] 

To  obtain  a  belter  understanding  of  chirality  and  assay 
its  future  role  in  electrical  design,  we  shall  examine  in 
the  following  pages  the  interaction  between  electro¬ 
magnetic  waves  and  chiral  objects  In  particular,  we 
shall  study  the  case  of  a  composite  medium  consisting 
of  randomly  oriented  chiral  conductors  embedded  in 
a  dielectric 


1.  Two  Conjectures  on  Chiral  Objects 

Chirality  is  a  purely  geometric  notion  which  refers  to 
the  lack  of  symmetry  of  an  object.  By  definition,  an 
object  is  chiral  if  it  cannot  be  brought  into  congruence 
with  its  mirror  image  by  translation  and  rotation  An 
object  that  is  not  chiral  is  said  to  be  achiral  Thus  all 
objects  are  either  chiral  or  achiral  Some  chiral  objects 
occur  naturally  in  two  versions  related  to  each  other  as 
a  chiral  object  and  its  mirror  image  Objects  so  related 
are  said  to  be  enantiomorphs  of  each  other. 

A  chiral  object  has  the  property  of  handedness :  it  must 
be  either  left-handed  or  right-handed  If  a  chiral  object 
is  left-(  right- fhanded.  its  enantiomorph  is  rtght-i  left  - 1 
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Fig  1  A  sketch  showing  chiral  objects  Heft  columnl  and  their 
enannomorphs  (right  columnl  From  top  to  bottom  are  thown  a 
helix,  a  Mobius  strip,  an  irregular  tetrahedron  and  a  glove 


handed  For  example,  if  the  chiral  object  is  a  left- 
( right- (handed  helix,  its  enantiomorph  is  a  nght-Ueft-l 
handed  helix. 

The  handedness  of  helices  was  made  clear  by 
Lindtnan's  and  Pickering's  experimental  results  which 
showed  that  a  collection  of  randomly  oriented  left- 
handed  helices  would  rotate  the  plane  of  polarization 
of  a  linearly  polarized  microwave  one  way  but  that  a 
collection  of  randomly  oriented  right-handed  helices 
would  rotate  the  plane  of  polarization  the  opposite 
way 

Assuming  that  this  relation  between  the  handedness  of 
the  helices  and  the  sense  of  rotation  of  the  microwave 
is  not  peculiar  only  to  helices  but  is  a  property  of  all 
chiral  objects  and  their  enantiomorphs,  we  are  led  to 
the  following  conjecture  Any  medium  composed  of 
randomly  oriented  equivalent  (simply-connected)  chi¬ 
ral  objects  will  rotate  the  plane  of  polarization  one 
way.  say.  to  the  left,  while  a  medium  composed  of  the 
enantiomorphs  of  these  objects  will  rotate  the  plane  of 
polarization  the  opposite  way,  i.e.,  to  the  right. 

In  Fig  1  we  see  common  examples  of  chiral  objects :  a 
helix,  a  Mobius  strip,  an  irregular  tetrahedron,  and  a 


2b  2b 

— H  (•—  — •<  H 


left-handed  right-handed 

Fig  2  Idealized  short  helices  used  in  calculations  The  plane  of  the 
loop  is  perpendicular  to  the  axis  of  the  straight  portion  of  the  «ire 


glove  On  one  side  of  the  Figure  is  the  chiral  object  and 
on  the  other  is  its  enantiomorph. 

A  type  of  (multiply-connected)  chiral  object  that  has 
recently  attracted  considerable  attention  is  the  wire 
braid.  The  theory  of  braids  is  a  developing  branch  of 
topology  [14.  15]  and  a  study  of  how  an  electromag¬ 
netic  wave  interacts  with  a  braid  may  help  in  the 
development  of  the  theory 

Examining  the  forces  that  are  exerted  on  certain  simple 
chiral  conFigurations  of  wire  when  an  electromagnetic 
wave  falls  on  them,  we  conjecture  that  the  forces  are 
such  as  to  reduce  the  chirality  of  the  configurations 
This  is  true  for  the  wire  helix,  for  the  three-stranded 
braid,  and  appears  to  be  true  in  general  This  tendency 
of  the  forces  makes  the  object  more  nearly  sym¬ 
metrical. 


2.  The  Short  Helix 

To  demonstrate  the  plausibility  of  the  above  conjec¬ 
tures.  we  examine  the  scattering  of  electromagnetic 
waves  from  a  metallic  chiral  object  For  computational 
simplicity,  the  chiral  object  is  chosen  to  be  an  electri¬ 
cally  small  perfect  conductor  having  the  form  of  a 
short  right-  or  left-handed  helix,  as  shown  in  Fig  2. 
The  calculation  is  simplified  by  referring  the  incident 
and  scattered  waves  to  the  scattering  plane  defined  by 
the  incident  and  scattered  wave  vectors  k  -ki,  and 
k  -  ki,,  respectively  (Fig  3).  The  incident  plane  wave  is 
composed  of  the  electric  field 

E‘ «■ +  a  e'V  le'u  III 

and  the  corresponding  magnetic  field 

B' *  E'  c. 


(2l 
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where  c  is  the  free-space  speed  of  light.  a  and  6  are 
real  numbers  with  a^+a:  =  1.  and  :  is  the  distance 
along  e„(=cj|xej.  The  circumflexed  quantities  are 
unit  vectors,  the  primes  denote  quantities  associated 
with  the  incident  wave,  and  the  subscripts  identify 
quantities  parallel  or  perpendicular  to  the  scattering 
plane.  The  harmonic  time  dependence  exp(-i<oM 
(where  u>  =  ck)  has  been  suppressed 
The  scattered  electric  Held  Ejffi  depends  on  the 
observation  angle  6.  defined  bs  the  relation 
cosfi  =  e.  e„.  and  on  the  induced  electric  and  magnetic 
dipole  moments  p  and  m  It  is  apparent  (Fig.  2]  that 
both  dipole  moments  are  directed  parallel  to  the  axis 
of  the  helix  This  axis  lies  along  the  unit  vector  e4 
whose  orientation  angles  are  i  and  /( (Fig  3 1 
The  incident  electric  field  induces  currents  in  the 
straight  portion  of  the  chiral  object,  and  bv  continuitv 
these  currents  must  also  flow  in  the  circular  portion  of 
the  object  The  current  in  the  straight  portion  contri¬ 
butes  to  the  electric  dipole  moment  of  the  object  and 
the  current  in  the  circular  portion  contributes  to  its 
magnetic  dipole  moment  In  a  complementarv  manner, 
the  incident  magnetic  field  induces  currents  in  the 
circular  portion  and  b>  continuitv  in  the  straight 
portions  Thus,  also  the  magnetic  field  contributes  to 
the  electric  and  magnetic  dipole  moments  of  the  object 
In  a  first-order  (Bonn  approximation  we  find  from  the 
heuristic  argument  above  that  the  electric  and  mag¬ 
netic  dipole  moments  of  the  ohiect  are  given  bv 

p  =  i„[/,.(.'J  E  I  ±  i/^ole,  B  l]<'j  (3l 

■*»/  B  i E  i ].„  (4i 

Flere.  as  in  (he  remainder  of  the  report,  the  upper 
llowerl  sign  corresponds  to  the  righl-handed  lleft- 
handedi  helix  of  Fig  2  The  permittivitv  the  per- 
meabtluv.  and  the  impedance  of  free  space  are  denoted 
bv/,,.  a,,,  tji  =(/i,,  i ,,  1 1  :l  The  electric  and  magnetic  self- 
susceptibilities.  /,  and  /„.  are  real  positive  quantities  as 
are  the  cross-susceptibilities  /rm  and  /mr  Clearlv.  / 
and  are  the  usual  electric  and  magnetic  susceptibil¬ 
ities  associated  with  electricullv  small  metallic  bodies 
I  he  cross-susceptibilities  /rm  and  are.  in  a  certain 
sense  a  measure  of  chiralitv  or  handedness  since  for 
achiral  bodies  -  0 

t  sing  known  approximation^.  the  self-suscepnhiUties 
can  he  written 


Fig  *  Vectors  indicating  (hr  directions  and  angles  of  incident  and 
scattered  waves  and  the  orientation  angles  of  the  chiral  object  The 
axes  are  oriented  in  such  a  manner  that  the  unit  vectors  » 

<*, .  «v  and  <„  are  in  the  scattering  plane  defined  b>  k  and  k 
The  unit  vectors  i  and  e  are  perpendicular  to  ihi^  plane 


arc  given  hv 

na2ki  |7| 

=  /,< 2/>  (Si 

From  phvsical  considerations  li  e  .freedom  from  ohm¬ 
ic  losses  and  conservation  of  energv  i.  n  appears  that 
/,„  and  are  equal  and  real,  i.c . 

/,m  =  /„.=/.  Ol 

where  /,  is  (heir  real  common  value 
li  follows  from  (5h9i  lhai  the  constraint 

/-(  =<■<  ;  dOi 

is  placed  upon  the  inductance  and  capacitance  of  the 
helix  and  the  common  value  /  for  the  cross- 

susceptibilities  is  related  to  the  inductance  and  capaci¬ 
tance  bv 

/  =  2/(nu;lt7  t  iL  =  2/(Jtu*l>;i  >C  1 1 1 1 


/,  =i2/r(  i„  (5i 

/„  =  irtw-'r'ji,.  /  (hi 

where  (  and  /  are  respectively  the  capacitance  and 
ihe  inductance  of  the  bodv  [lb],  and  2/  and  2u 
represent  the  length  and  the  width  of  the  short  helix 
1 1  ig  2 1  It  can  be  shown  that  the  cross-susceptibilities 


From  the  knowledge  of  p  and  m  the  scattered  field  can 
be  calculated  bv  the  formula 


„  k 

E,  (ifi=  [((■«■(»('-<'«  m  <  1  ( 1 2 1 

4m,,r  ”  ” 

To  gain  further  insight  into  the  problem,  it  o  useful  t >> 
find  the  constitutive  relations  of  a  medium  composed 


214 


D  L  .laggard  el  al 


of  randomly  oriented  equivalent  chiral  objects.  These 
consitutive  relations  must  have  the  form  [17] 


the  scattered  field  is  considered. 


P=,E  +  „.,B 
M  = '. -,E  +  y„B . 


(13) 

(14) 


lda(9)\ 

\dQ  , 


(k2S'J\)2  , 

-2 —  IZ,  —  1  -vcos0)- 


1024rr 


(21) 


where  P  and  M  are,  respectively,  the  polarization  and 
magnetization  of  the  medium. 

Energy  conservation  dictates  that  for  a  lossless 


when  the  incident  wave  is  left  circularly  polarized 
(LCP)  and  only  the  left  circularly  polarized  part  of  the 
scattered  field  is  considered,  and 


medium 

(d<Ji9)\ 

/dm\ 

_ _ .* 

1  me  >  rm  ' 

(15) 

'  ^2  /RCP  — LCP 

V  dQ  1  LCP-*CP 

where  the  asterisk  denotes  complex  conjugate.  If 
and  ,rm  not  only  satisfy  (15)  but  also  are  purely 
imaginary  quantities,  then  the  constitutive  relations 
(13)  and  (14)  are  those  of  an  optically  active  medium 

cm 

To  find  the  consitutive  parameters  for  a  medium 
composed  of  N  non-interacting  short  helices  per  unit 
volume,  we  compare  (3)  and  (4),  averaged  over  orien¬ 
tation  angles  x  and  p,  with  (13)  and  (14).  Thus  we 
obtain 


V.=  Vw.,/4 

7rm= 

:mr=+iNxmr,*’l 


(16) 

(17) 

(18) 
(19) 


Since  x„  and  are  real  and  equal,  we  see  from  (18) 
and  (19)  that  •/„,  and  ytm  are  purely  imaginary  and 
satisfy  (15).  Hence  the  medium  composed  of  short 
helices  exhibits  optical  activity. 


(*2NdK) 


10247T5-  lz'  +  *"1' (1  "cose»  <22> 

when  the  incident  wave  is  RCP  or  LCP  and  the 
scattered  field  is  LCP  or  RCP.  respectively. 

From  these  expressions  we  see  that  for  scattering  in  the 
forward  direction  (6  =  0)  right  (left)  circular  polar¬ 
ization  produces  a  right  (left)  circularly  polarized  scat¬ 
tered  field  whereas  in  the  backward  direction  (6  =  it) 
right  (left)  circular  polarization  produces  a  left  (right) 
circularly  polarized  scattered  field 
Next,  we  consider  a  wave  normally  incident  on  an 
electrically  thin  slab  of  width  d  which  contains  S' 
randomly  oriented  helices  per  unit  volume  Using 
again  the  assumption  that  the  helices  are  non- 
interacting,  we  find  from  averaging  (12)  that  the  trans¬ 
mitted  field  is  given  by 


E,r  = 


Vn  +  fl-c  v- 


iNkdL  ,  4_. 


3.  Collection  of  Short  Helices 


+  i[±Vxf-ifl-ci*(X*-X»)^- 


To  find  the  scattered  field  of  a  collection  of  randomly 
oriented  identical  helices  we  can  use  one  of  two 
approaches.  One  approach  uses  (12)  averaged  over 
orientation  angles  x  and  P:  the  other  approach  uses 
(13l  and  (14)  directly.  These  two  approaches  give  the 
same  result.  Here  we  use  the  former  of  the  two 
approaches 

Let  us  suppose  that  we  have  a  collection  of  N  non- 
interacting  helices  per  unit  volume  occupying  a  small 
volume  AV  When  the  incident  wave  is  circularly 
polarized,  the  scattering  cross-sections  per  unit  solid 
angle  Q  are  found  to  be 


/d( 7(0) 


PCP  -PCP 


(k2NJV)J 

1024r: 


Iz,-Xw±2xfl2(l  +cos  0): 


(20) 


when  the  incident  wave  is  right  circularly  polarized 
i  RCP)  and  only  the  right  circularly  polarized  part  of 


and  the  reflected  field  by 

Er,r-  ^y^(Z,  +  X»»[0||f|i  +  ‘J.e,'e  ]e',t:  (24) 

where  the  above  expressions  are  correct  to  the  first 
order  in  ( Nkdx )  (here  x  stands  for  Xm  or  X,l 
From  (23)  it  can  be  shown  that  the  plane  of  polar¬ 
ization  is  rotated  through  the  angle  <t>  where 

<t>  s  tan  <t>  =  T  NkdX'/d  (25) 

for  waves  which  pass  through  the  chiral  medium  Here 
<t>  is  measured  from  eH  towards  e  Expression  (25)  is 
again  correct  to  the  first  order  in  (A’Ady).  This  equation 
expresses  a  genera!  result  which  holds  for  any  medium 
composed  of  objects  characterized  by  parallel  electnc 
and  magnetic  dipole  moments  with  non-zero  cross¬ 
susceptibilities  For  the  short  helices  pictured  in  Fig  2 
we  can  find  a  lower  bound  on  the  capacitance  C  by  the 
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Fig  4  A  schematic  of  a  short  helix  (chiral  1  evolving  into  a 
straight  wire  (achiral  I  under  the  influence  of  forces  produced  h> 
Li  induced  currents 


expression  [18] 

CSfc0(47t):  3  (26) 

where  lk(  =  2nh2(/  +  rw)l  is  the  volume  occupied  by  the 
short  wire  helix  and  h  is  the  wire  radius.  This  in¬ 
equality.  with  the  aid  of  (11)  and  (25|.  yields  the 
following  lower  bound  for  the  magnitude  of  the  ro¬ 
tation  angle 

|d>|  i  |tan  01  §4' 1  iit2  3N(MH2/Mmj2X3fcH'(1l1  V  (27) 


g  fS 

Ls> 


This  bound  is  proportional  to  the  product  of  the  third 
root  of  the  volume  fj,  of  the  wire  helix  and  the 
cylindrical  volume  containing  the  helix  (=2 Ina2). 
From  (23)  the  eccentricity  of  the  transmitted  polar¬ 
ization  ellipse  differs  from  that  of  the  incident  polar¬ 
ization  ellipse  by  a  factor  of  order  I Nkd /)2  However, 
this  transmitted  field  is  correct  only  to  order  (A 'kd/l 
Therefore,  the  change  in  the  eccentricity  of  the  polar¬ 
ization  ellipse  cannot  be  determined  exactly  from  this 
model.  To  the  first  order  in  (A kd/)  the  eccentricity  is 
unchanged 

The  reflected  wave  (24)  to  order  {Nkd/)  shows  zero 
rotation  for  the  plane  of  polarization  and  zero  change 
in  the  eccentricity  for  the  polarization  ellipse 
Therefore,  for  reflected  waves,  the  slab  of  chiral  me¬ 
dium  behaves  as  an  ordinary  dielectric  slab.  These 
polarization  characteristics  are  due  to  the  fact  that  in 
the  backscatter  direction,  in  the  first  order,  the  effects 
of  chirality  are  not  present  in  the  scattered  field  [see 
(20).  1 2 1 1.  and  (23 1] 

From  Noether's  theorem  [19]  it  can  be  shown  that  the 
angular  momentum  of  the  electromagnetic  field  is 
conserved  for  a  medium  described  by  (13H19)  This 
implies  that  no  torque  is  exerted  on  a  slab  of  chiral 
medium  It  is  not  surprising  that  there  is  no  torque 
since  the  electrical  properties  of  the  slab  are  invariant 
under  rotations  of  the  slab  about  (•„  With  a  knowledge 
of  the  state  of  polarization  of  the  incident  wave, 
conservation  of  field  angular  momentum  further  im¬ 
plies  that  the  state  of  polarization  of  the  reflected  wave 
can  be  determined  from  the  state  of  polarization  of  the 
transmitted  wave,  and  vice  versa  Some  experimental 


Fig  5  Under  ihe  action  of  forces  produced  b>  the  induced  currents  a 
long  helix  (shown  at  left)  is  graduallv  shortened  into  more  closelv 
spaced  loops  of  increased  radius  (shown  at  righn 

results  indicate  that  there  is  a  change  of  eccentricity 
between  the  incident  and  scattered  fields  due  to  the 
chiral  medium  [II.  12]. 

From  the  above  considerations,  the  conjecture  that  a 
collection  of  chiral  objects  will  rotate  the  plane  of 
polarization  becomes  plausible 

4.  Reduction  of  Chirality 

Assuming  that  the  helix  in  Fig  4  is  made  of  flexible 
wire,  we  see  that  the  currents  that  are  induced  tend  to 
deform  the  helix.  The  current  along  the  circular  por¬ 
tion  tends  to  open  up  the  circle  and  make  a  planar 
figure  out  of  the  original  helix  Moreover,  interaction 
with  the  current  along  the  straight  portion  of  the  helix 
tends  to  elongate  the  planar  figure  into  a  straight  line 
Since  planar  figures  are  achiral,  we  thus  see  that  the 
helix  evolves  into  a  planar  figure  and  that  the  chirality 
of  the  configuration  is  reduced 

Suppose  now  that  we  have  a  flexible  helix  of  manv 
turns  (Fig  5l  In  this  case  the  induced  current  forces 
adjacent  turns  together  and  at  thv  same  time  makes 
each  turn  expand  into  a  turn  of  larger  radius  Thus  the 
original  helix  becomes  a  shortened  heltx  of  larger 
radius  Since  the  shortened  helix  is  less  chiral  than  Un¬ 
original  helix,  we  see  that  here  again  the  induced 
currents  tend  to  reduce  the  chirality  of  the  con 
figuration 
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Fig  6  Bus  bars  and  L:  connected  b>  three  strands  of  the  braid 
defined  b>  n :nt  '  a  > 


Another  type  of  chiral  object  is  a  braid  of  non- 
intersecting  wires  Following  Artin's  theory  of  braids 
[  1 4]  we  may  describe  a  braid  by  projecting  it  on  a 
plane  and  expressing  the  projected  pattern  as  the 
product  of  terms,  each  of  which  is  n,  or  a~  1  Here  a, 
denotes  that  the  strand  in  position  i  crosses  in  front  of 
the  strand  in  position  i  +  1  and  it  1  denotes  that  the 
latter  crosses  in  front  of  the  former 
Let  us  consider  the  three-stranded  braid  <T,<rj  'o: 
shown  in  Fig  6  The  bus  bars  and  L.  are  connected 
by  three  flexible  wires  at  the  freely  movable  but 
ordered  terminals  Pr  P:.  P,  and  Qr  Q:.  Q ,  Clearly, 
the  braid  and  bus  bars  form  a  chiral  object  An 
incident  wave  will  induce  currents  in  the  braid  and  bus 
bars  and  these  currents  will  deform  the  configuration, 
wz..  will  make  it  more  nearly  planar,  and  thus  reduce 
its  chirality 

From  these  and  other  special  cases  we  are  led  to  the 
conjecture  that  the  electrodynamic  forces  are  such  as 
to  reduce  or  minimize  chirality 

Conclusions 

\ke  direct  attention  to  the  interaction  of  electromag¬ 
netic  fields  with  macroscopic  chiral  objects.  By 
examining  the  wire  helix  and  the  wire  braid  as  chiral 
objects,  we  obtain  results  which  conform  to  the  conjec¬ 


ture  that  composite  media  composed  of  macroscopic 
chiral  objects  are  optically  active,  and  to  the  conjecture 
that  electrodynamic  forces  tend  to  reduce  chirality 
These  considerations  are  expected  to  play  a  role  in  the 
development  of  diagnostic  tools  for  remote  sensing,  in 
the  design  of  electromagnetic  shields  and  in  the  pre¬ 
diction  of  structural  deformations. 
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Steady  State  and  Transient  Electromagnetic  Coupling  Through  Slabs 
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Abstract— The  problem  of  electromagnetic  transmission  through  a 
slab  where  transmitting  and  receiving  antennas  are  at  finite  distances 
from  the  slab  is  considered.  The  mathematical  formulation  of  the 
problem  is  quite  general.  A  detailed  solution  is  presented  for  the 
case  of  a  highly  conducting  slab  exposed  to  sinusoidal  and  transient 
excitations.  A  discussion  is  given  of  the  conditions  under  which  meas¬ 
urements  with  the  source  and  receiver  at  finite  distances  are  equivalent 
to  ihe  same  measurements  with  plane  wave  excitation. 

1  INTRODUCTION 

NE  OF  THE  simplest  conceivable  ways  for  determing  the 
electromagnetic  properties  of  materials  is  to  measure  the 
electromagnetic  field  transmitted  through  a  slab  of  the  material 
under  test.  The  corresponding  mathematical  model  consists 
of  an  infinite  slab  with  transmitting  and  receiving  antennas 
placed  on  opposite  sides  of  the  slab.  The  model  provides  a 
reasonably  good  approximation  to  the  real  situation  of  a  slab 
of  finite  extent  when  the  distance  between  transmitting  and 
receiving  points  is  small  compared  to  the  transverse  slab 
dimensions. 

Measurements  can  be  made  in  the  sinusoidal  or  the  tran¬ 
sient  regime.  For  instance,  MIL  standards  for  evaluating  the 
shielding  effectiveness  of  materials  [  1  ]  require  that  transmis¬ 
sion  measurements  be  made  in  the  steady  state  at  prescribed 
frequencies  and  then  in  a  pulsed  regime  using  wire  and  loop 
antennas  placed  at  prescribed  distances  from  the  slab  of  shield¬ 
ing  material  Although  these  standards  are  useful  for  relative 
comparisons,  a  fundamental  question  remains  unanswered 
does  the  measurement  depend  only  on  the  electromagnetic 
properties  of  the  slab  (and  on  its  thickness),  or  does  it  depend 
also  on  antenna  type  and  orientation,  antenna  distance,  and 
(for  transient  measurements)  on  transmitted  waveform? 

A  crude  but  simple  method  for  studying  (or,  at  least,  having 
an  estimate  of)  the  field  coupled  to  the  inside  of  an  enclosure 
is  to  consider  the  transmission  through  a  slab,  provided  the  en¬ 
closure  is  large  in  terms  of  the  incident  wavelength  The  slab 
ma>  be  perforated,  inhomogeneous,  or  described  by  stochastic 
parameters,  the  last  case  being  relevant  to  near-millimeter  prop¬ 
agation  through  aerosols  used  for  camouflage  tactics.  In  electro¬ 
magnetic  pulse  (EMP)  experiments  it  is  customary  to  simulate 
the  EMP  plane  wave  signal  by  using  rather  sophisticated  an¬ 
tennas  and  guiding  devices  ( 2) ,  {31 .  An  attractive  alternative 
to  this  approach  can  result  from  an  understanding  and  ex¬ 
ploration  of  the  role  played  by  localized  sources  at  finite  dis¬ 
tances  from  the  test  object 
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The  objectives  of  this  paper  are  to  reconsider  the  problem 
of  steady-state  and  transient  coupling  through  a  slab  with 
transmitting  and  receiving  antennas  located  at  finite  distances 
from  the  slab,  to  cast  the  problem  in  an  elegant  form;  and  to 
show  that,  at  least  in  the  case  of  a  highly  conductive  slab,  sim¬ 
ple  analytical  solutions  to  the  problem  can  be  obtained  An 
important  result  of  the  paper  is  the  determination  of  antenna 
positions  and  (in  the  tiansient  regime)  of  incident  waveforms 
that  will  yield  a  transmitted  field  practically  the  same  as  that 
produced  by  plane  wave  excitation. 

Transmission  through  highly  conductive  slabs  is  certainly 
not  a  new  problem.  For  plane  wave  steady-state  excitation, 
transmission  line  techniques  can  easily  be  applied  |4)  For 
pulsed  plane  wave  excitations,  the  solution  is  also  available  [5] 
The  situation  is  much  less  satisfactory  for  the  case  we  want  to 
study.  It  is  not  the  aim  of  this  paper  to  provide  a  full  biblio¬ 
graphy  on  this  subject  (for  a  more  complete  bibliography  see 
[6) ).  We  note  only  that  the  first  attempt  to  solve  this  problem 
was  made  in  1936  1 7 )  by  accommodating  the  classical  results 
of  Maxwell  on  eddy  currents  and  thin  shields  to  the  case  of 
two  coaxial  loops  separated  by  a  plane  conducting  sheet  Early 
studies  on  antenna  coupling  through  plane  shields  were  based 
on  low-frequency  18],  |9]  or  quasi-static  [  10]  approxima¬ 
tions.  were  mainly  relative  to  loop  excitation  |8)-|)0).  and 
required  numerical  computation  [  8 )  — [  12}  of  integral  expres¬ 
sions  for  the  transmitted  field.  Although  the  validity  of  the 
simple  transmission  line  theory  [4)  for  antennas  at  finite  dis¬ 
tances  from  the  shield,  or  shields  of  finite  extent,  has  been 
questioned  |  13),  it  appears  that  all  expressions  derived  in  the 
referenced  literature  resemble  Schelkunoffs  formulas  I  14] 

Due  to  the  symmetry  of  the  problem  it  can  easily  be  sur¬ 
mised  that  plane  wave  expansion  techniques  provide  a  power¬ 
ful  tool  of  analysis  for  an  arbitrary  type  of  excitation  o:  ar, 
infinite  slab.  These  techniques  have  been  recently  applied  1 14) 
[15]  to  the  case  of  electric  or  magnetic  dipole  excitation  in 
parallel  (dipoles  parallel  to  the  slab)  or  coaxial  (dipoles  normal 
to  the  slab)  configuration,  by  computing  the  transmitted  field 
through  the  use  of  fast  Fourier  numerical  programs  In  this 
paper  we  shall  use  the  same  approach  However,  we  will  show 
that,  although  the  Fourier  transformation  of  the  fields  is  a 
logical  intermediate  step  of  the  analysis,  it  is  not  needed  in  the 
final  formulation  of  the  solution  Indeed,  the  solution  can  be 
conveniently  expressed  in  terms  of  a  convolution  integral, 
wherein  the  presence  of  the  slab  is  described  by  an  appropriate 
transfer  function.  Then,  at  least  for  antennas  in  coaxial  con¬ 
figuration.  the  convolution  integral  can  be  analytically  evalu¬ 
ated  both  in  steady-state  and  transient  regimes,  and  no  numer¬ 
ical  work  is  necessary.  Inspection  of  the  solution  allows  us  to 
answer  the  original  question  about  the  influence  of  the  finite 
antenna  separation  on  measurements.  After  all  the  mathemat¬ 
ical  machinery  has  been  worked  out  and  simple,  physically 
sound,  understandable  results  are  obtained,  a  discussion  of  the 
final  results  is  presented  in  Section  VI 
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11  CIRCUIT-LIKE  ANA?  VS1S  OF  ELECTROMAGNETIC 
TRANSMISSION  THROUGH  A  SLAB 

With  reference  to  Fig.  1.  let  us  consider  an  infinite  slab  of 
thickness  s  and  characterized,  in  the  frequency  domain,  by 
permittivity  e  =  e0e,.  permeability  /t  =  u 0nr,  and  conducti¬ 
vity  o  We  want  to  compute  the  field  E(,  H*  transmitted  at 
:  >  s  along  the  :  axis  when  the  incident  field  E1,  H\  i.e.,  the 
field  produced  by  the  sources  when  the  slab  is  removed,  is 
known  at  :  =  0.  For  this  purpose  it  is  convenient  to  expand 
the  incident  field  in  a  plane  wave  set.  since  the  mteraction  of 
individual  plane  wave  components  with  the  slab  can  be  con¬ 
veniently  taken  into  account. 

Accordingly,  let  H,‘ f.r,  y,  01.  t,‘ f.r,  y.  0)  be  the  r  com¬ 
ponents  of  the  field  incident  on  the  slab  surface,  with  an  as¬ 
sumed  time  dependence  exp(/cor)  The  corresponding  spectra) 
components  ht‘(u,  v),  e,'(u,  v)  are  given,  at  z  =  0,  by 

h,‘(u.  v)  =  — l—  [  dx  (  dy  y.  0) 

<2rr)J 


wherein  x  =  ux  +  vy  +  uf  and  is  the  propagation  vector  re¬ 
ferring  to  a  Cartesian  system  of  unit  vectors  x,  r,  and  the 
upper  (lower)  signs  refer  to  waves  propagating  in  the  positive 
(negative)  direction  of  the  z  axis.  Equations  (2.S)  and  (2.6) 
represent  the  total  spectral  field  as  a  superposition  of  TE  (e,  = 
0)  and  TM  (h,  =  0)  parts.  With  the  medium  being  identical  at 
both  sides  of  the  slab,  it  is  then  evident  that  r»  coincides  with 
the  usual  slab  transmission  coefficient  for  TE  plane  wave  in¬ 
cidence  and.  similarly,  iE  is  the  same  as  the  slab  transmission 
coefficient  for  TM  plane  wave  incidence.  Letting 


W  =  y/K2  —  ( u 1  +  P2)  , 


urw  a  +  /u£0er  w 


we  have 


r(u.  v)  = 


exp  (-/wfr) 


<  i  +  y)2 


i  - 


i  -tV 


i  +  y 


(2.9) 


exp  (-2;w,r) 


wherein  7  may  be  taken  equal  to  yH  or  yE  in  order  to  obtain 
1 H  or  r£.  respectively,  and  x  =  un/foPo 

The  spectral  components  hJt  ez  at  any  r  >  r  are  equal  to 
the  corresponding  values  (2.3)  and  (2.4)  at  r  =  1  tunes  the 
plane  wave  transfer  function  exp[ -jw(z  -  r)) .  Accordingly, 
the  r  component  Ft‘ (x,  y,  r)  of  the  field  transmitted  at  any 
arbitrary  abscissa  :  >  j  will  be  expressed  in  terms  of  the 
double  Fourier  integral 


•  exp  (/ux  +  /vy) 


e,'(u.  V) 


- -  [  dx  f  dy  £V(x.  y.  0) 

(2ff)J  J-m  J-. 


(2.1  I 


Ft  '(x.  y,  r)  = 


dvft'(u,  ir)r(u,  v) 


•  exp  [-/w(i-r))  *  exp  (-/ux -/vy)  (2.10) 


•  exp  (/ux  +  ivy)  (2.2) 

At  :  =  r.  i.e  .  at  the  output  of  the  slab,  the  spectral  com¬ 
ponents  h/(u,  it).  f,'(u.  v )  will  be  linearly  related  to  the  inci¬ 
dent  components  (2.1)  and  (2.2)  in  the  case  of  a  slab  made  of 
a  linear  material  Hence, 

h/(u.  v)  =  iH(u.  V)hj'(u.  v)  (2.3) 

e/(u.  1;)  «  tE(u.  v)t,‘(u.  V)  (2.4) 

The  transfer  coefficients  tH.  rE  can  be  easily  computed  for  a 
homogeneous  isotropic  slab  by  noting  that  the  transverse 
spectral  components  h,(u,  it),  t,(u,  it)  are  related  to  the  longi¬ 
tudinal  ones  h,(u,  v).  c,(u.  u)  via  the  following  relations: 


h,= 


wee,xx2  ♦  wh,:x(Kxi) 
u2  +  v2 


(2.5) 


CJJlhjKxi  *  Wft2x(xxi) 

U1  +  It* 


(2  6) 


wherein  may  be  taken  equal  to  h,‘  or  e,‘  and,  correspond¬ 
ingly,  the  values  of  iH  or  tE  should  be  used. 

On  the  other  hand,  the  spectral  representation  of  the  r 
components  of  the  incident  field  (the  slab  is  now  removed)  at 
any  abscissa  c  is  obviously  the  following: 


F,‘(x.  y.  z)  — 


dv  fi‘(u.  v)  exp  (— fwz) 


•  exp(-/ux -/try).  (2.11) 

Comparison  of  (2.10)  and  (2.11)  shows  that  the  transmitted 
field  can  be  computed  as  the  double  convolution  of  the 
incident  field  and  the  double  Fourier  transform  of  r(u,  v) 
exp (jws).  hence 

F,'(x.  y.  r)  =  J  dx' J  dy'  F,‘(x‘.y',  r) 

•  T(x-x'.y-y')  (2  12) 
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dl'  r(u,  v )  exp  (;ns) 


•  exp  (—  jux  —  ivy) 


(2.13) 


In  the  words  of  system  theory'  F1  is  identified  with  the  output 
of  a  linear  system  described  by  the  unit  response  function 
(2.13)  and  excited  by  the  input  F‘. 

We  further  note  that  relations  similar  to  (2.12)  exist  be¬ 
tween  the  transmitted  and  incident  transverse  components  of 
the  field,  as  easily  follows  from  (2.5)  and  (2.6).  It  is  only  ne¬ 
cessary  to  decompose  the  incident  field  in  its  TE  and  TM  parts 
and  then  to  apply  superposition. 

III.  THE  AZIMUTH  ALLY  SYMMETRIC  CASE 

A  case  of  particular  interest  is  obtained  when  the  incident 
field  does  not  depend  on  x  and  y  separately  but  rather  upon 
the  transverse  coordinate  p  =  \Jx2  +  y 2  For  instance,  if  the 
source  is  taken  equal  to  an  elementary  electric  or  magnetic 
dipole  parallel  to  the  :  axis  at  PIO, 0,  -d),  then 

F,'(x.y.  z)  =  F,'(p.  .-)  =  --^-[k2A+  77  •  A]  •  f 

K.1 


A 


\J 

i 

* 


Fig  2.  Integration  path  in  complex  u  plane. 


wherein  the  order  of  integration  has  been  reversed.  Upon  sub¬ 
stitution  of  (3.1)  in  (3.51  the  inner  integral  can  be  evaluated 
by  repeated  integration  by  parts  as  follows: 


i: 


P  dp  F‘(p.  d  +  2)J0(Pi ) 


;uC  ,  f  P  exp  ( ~)k\/p 2  +  (d  +  r)2) 

~(2  - ;  n=  . ft - Jo(Pi)dp 

Jo 


Vp2  +  (d  +  r)1 


i  OJ 


+ 


d2A~ 

d:2 


ICO  1  d  BA 

— - P—  (3.1) 

k2  p  dp  dp 


co£  exp  [-/n/k*  ~  ^(d  +  J)] 

k2  \Jk i  —  £* 


wherein 


A(p.:)  = 


^  exp  (— /kVp2  +  (d  +  2)2) 
\'pi  +  (d  +  2)^ 


(3.2) 


is  an  electric  or  magnetic  vector  potential,  the  source  intensity 
being  proportional  to  the  constant  C. 

The  integrals  (2.12)  and  (2.13)  can  now  be  simplified  by 
using  the  change  of  coordinates: 


the  last  expression  stemming  from  a  known  Fourier-Bessel 
transform  [  16) .  Note  that  \/k2  -  £2  =  -rJ £2  -  for  £2  > 
k2  and  that  we  have  implicitly  assumed  k  #  0  in  this  section. 

The  formal  expression  for  the  r  component  of  the  field 
transmitted  through  the  slab  is  now  the  following: 

•  F.HQ.  O.c)  =  wCk  /  ( 1  —  u2)r(u)  exp  (— /x/u)  du  (3.7) 
JT 


x=p  cos<>.  >'  =  psind),  u  =  £  cos  i^,  b=£sinv. 
Accordingly. 

Tlx.  y)  =  T(p)  =  /  {  d|  r(£ )  exp  (/V*2  -  £2r) 


(3.3) 


■i; 


i: 


exp  [  jp  \  cos  ( ip  <P)J  dy 

=  >/  $  d| /0(pi)'(*  >  exp  (/V*1  -  E4!).  (3.4) 

•'o 

and  the  field  transmitted  on  the  axis  is  given  by 

FZU 0,0.:)=  I  p  dp  Fxl(p.  d  +  z)T(p) 

'o 

=  f  £  d(  exp  (/V*2  - 
•'o 


f 

Jo 


p  dp  Fj'(p.  d  +  2)J0(ip)  (3.5) 


wherein  /  =  d  -  i  +  z,  the  integration  path  T  is  depicted  in 
Fig.  2,  and  the  substitution  k2  -  f2  =  k2u2  has  been  used. 

IV.  THE  CASE  OF  AN  ELECTRIC  PLANE  SHIELD- 
STEADY-STATE  EXCITATION 

A  case  particularly  interesting  for  applications  is  obtained 
when  pr  =  I,  c  >  co£0er,  i.e.,  when  a  highly  conducting  non¬ 
magnetic  slab  is  used  as  a  shielding  screen.  As  already  noted  in 
Section  I,  this  is  an  important  configuration  in  shielding  theory 
and  practice.  The  solution  to  this  problem  js  available  in  nu¬ 
merical  form  |  14) ,  [  15)  for  prescribed  sinusoidal  time  varia¬ 
tion  and  arbitrary  spatial  dependence  for  the  fields,  and  in 
analytical  form  [5)  for  prescribed  plane  wave  excitation  and 
arbitrary  time  variation. 

The  case  of  a  magnetic  dipole  excitation  is  considered  first. 
The  expression  for  r(u)  pertinent  to  this  case  is  the  following: 

4uV®2  +  u2 

'"(U)  =  (u  +  Vo^+u1)2 


exp  (— /v/ct^  +  u2kj) 


f  u  —  v’a2  +  u2  1 

2 

1  - 

_U  +  yja2  +  u2 

exp  (—2 is/a2  +  u^ks) 

(4.1) 


-set 
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„  a  ;cj«0(e,-  1  >  o 

a2  = - s -  <42) 

/  uie0  /we  o 

It  is  noted  that  rH(u)  exhibits  no  singularity  in  the  lower  right 
quadrant  of  the  complex  u  plane,  so  that  the  integration  path 
P  can  be  freely  deformed  therein,  e  g.,  in  the  new  path  f  (see 
Fig  2)  When  (4.1)  is  substituted  in  (3.1)  it  is  noted  that  we 
can  neglect  u2  with  respect  to  a2  provided  the  integrand  is 
negligible  when  u  >  la  I.  Accordingly,  when  k1  I  a  l>  I  the  in¬ 
tegral  (33),  specified  to  the  case  at  hand,  becomes  ■ 


4  exp  (— i<xks) 

H,'( 0.  0. 1)  =  — /wCx - exp  (— /kI) 

a  1  —  exp  (— 2/axj) 


v(—  /It2  +  3v  +  2/)  exp  (— kIv)  du  (4.3) 


and  the  origin  of  coordinates  is  now  in  correspondence  to  the 
source.  The  integral  is  now  straightforward  to  evaluate  and  can 
be  conveniently  normalized  to  the  value  of  the  incident  field 
H,* (0,  0, 1).  We  have 


H,‘( 0.  o.l) 


(” 4  exp  (—/oms) 


3 

1  +  —  + 

IK  l 


3 


H,‘( 0.  0.  /) 


a  1  —  exp  (— 2/o«a)J 


1 

1  +  — 

;«/ 


=  r0(a.  K!)nH(K() 


(4,4) 


It  is  noted  th*t  the  first  bracketed  term  r0(a,  xj)  is  just  the 
plane  wave  transmission  coefficient  under  normal  incidence 
and  appropriate  to  a  highly  conducting  screen.  The  second 
term  fi(x/)  depends  on  the  mutual  distance  /  between  trans¬ 
mitting  and  receiving  points  and  approaches  1  when  .£/  >  1. 
Accordingly,  it  follows  that  a  simple  plane  wave  transmission 
coefficient  can  be  used  for  evaluating  shielding  effectiveness 
provided  that  transmitting  and  receiving  antennas  are  a  few 
wavelengths  apart 

On  the  other  hand,  when  kI  is  small.  nH(x/)  s  3//k(.  and 


//,'( 0,  0. /) 

H,‘ '(0.  0. 1) 


3 

a  —  r0(o.  xj)  — 

IK  I 


exp  ( — /cxkj  —  /tt/4)36 
( 1  —  exp  (— 2/oxj)]\/T/ 


(4.5) 


wherein  6  is  the  skin  depth  of  the  screen.  Note  that  <4.5 )  is 
valid  provided  that  6/1  <  1 ,  otherwise  the  assumption  x  I  a  1/  > 
1  is  no  longer  met 

The  case  of  an  electric  dipole  excitation  can  be  treated  sim¬ 
ilarly  We  have 


4a2u\/a2  +  u2 
(a2u  +  u  +  Vcr5  +  u 2)2 

exp  (—/s/a*  +  u2ks) 


exp  (— 2/y/a2  +  u4xi) 

(4.6) 


t  _  I °2u  4  u  ~  y/a2  Ftfl  | 
\  a2  u  +  u  +  v  a2  +  u2  / 


We  can  now  neglect  u2  with  respect  too2  without  serious  lim¬ 
itation  in  the  validity  of  the  results.  The  integral  corresponding 


to  (4.3)  is  the  following: 

I  f  exp  (-/xfu) 

£,'(0.  0,  /)  =  cjCxt0(o.  ks)  •  j  /  -  du 

|  A  “ 


+  /J  1 1  -  /u)  exp  (— Klv)  dV  j  (4  7) 


which  can  be  easily  evaluated  to  yield 

£,'(0.  0. /) 

- =  t0(o.  Ks)ne(Kl) 

E,‘(0.  0. 1) 

1 +/x/  +  (jx/)2exp(jx/)(Ci(x/)-/Si(x/)] 

n£(x /)  =  $ - (4  s > 

i  +  — 

IK  I 

wherein  the  cosinus  integral  Cifrj  and  sinus  integral  Si(x) 
functions  [17]  do  appear. 

It  is  again  noted  that  fl£(x/)  -»  1  when  Kl  >  1,  as  easily 
follows  upon  use  of  the  asymtotic  series  expansions  [  1^]  of 
the  functions  Ci(x)  andSi(x),  so  that  (4.8)  reduces  again  to  the 
plane  wave  transmission  coefficient  i0(a,  ks)  provided  that 
transmitting  and  receiving  antennas  are  a  few  wavelengths 
apart  On  the  contrary,  when  x/  is  small,  a  proper  series  ex¬ 
pansion  [  17)  shows  that  fl£(  kI)  5s  yx(/2  and 

£,'(0.0./)  ik  l 
- a  —  :0(a.  ks) 


exp  (-)txxi  +/7T/4)  6 

=  / - (x/)2  —=^.  (4.9) 

1  —  exp  (— 2/axj)  s/ll 


V  THE  CASE  OF  AN  ELECTRIC  PLANE  SHIELD- 
TRANSIENT  EXCITATION 

We  have  shown  in  Section  IV  that  the  steady-state  z 
components  of  the  field  transmitted  through  a  highly  conduc¬ 
tive  plane  shield  are  given  by 

£,'(0.  0.  /)  =  £,'(0.  0.  /)r0(o.  xj)J2(x/).  (5  1) 

It  is  then  evident  that  the  r  component  of  the  transient  trans¬ 
mitted  field  can  be  obtained  by  time-convolving  the  transient 
z  components  of  the  incident  field  with  the  inverse  Fourier 
transforms  of  t0(u>)  and  fi(u>),  let  us  say  T0(t)  and  H(r). 
Use  of  Laplace  inversion  tables  [18]  shows  that 


exp  (— n2r;/r) 

,5/2 


(2n27)  —  r) 


£_  exp  (— n2T)/t) 
dt  i */* 


d 

~  so (0 
dt 


(5.2) 


where  r  =  e0lo  and  is  the  relaxation  time  of  the  material  of 
the  shield,  and  tj  *  j2/(c2  r)  and  is  the  diffusion  time  through 
the  shield  thickness. 


4  j 
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F  ig  3  Qualitative  behavior  of  first  series  lerm  ol  lunclion  TqU) 


I  it  4  Qualitative  behavior  of  pulsed  field  ol  time  duration  T  alter 
transmission  through  a  highly  conductive  slab 


A  qualitative  behavior  of  the  first  term  n  =  1  of  70tn  is 
given  in  Fig  3,  wherein  =  4  \Z~!2iU'r}  =  S 0<  3tji  and  is 

the  maximum  value  of  the  function  50(rl  (see  also  Fig  4). 
where  in  ,5mi,  "e"  is  the  Neper's  constant 

The  behavior  of  successive  terms  of  the  series  (5 .2  I  is  similar 
to  that  depicted  in  Fig  3  The  maxima  occur  at  later  times  and 
their  absolute  values  arc  smaller  by  the  factor  expl  -  3  61  n- 
1  l]  n Accordingly  .  they  can  be  safely  neglected,  and  we  can 
taise  only  the  first  term  ol  the  series  (5  3  1 . 

Afier  some  algebra  Laplace  inversion  [HI  of  the  two  func¬ 
tions  5’l oc I  leads  to 

3  -  exp  (-  r  T ) 

r. bd)  + - IV}  i  5  3  i 


a£uj  =  in i - 


T 


exp  I 


:'T)LV) 


exr  (-1/) 


i54i 


where  6(rl  and  (  (fl  are  the  Diras  and  the  unit  step  lunction. 
respectively,  T  -  .  <  and  is  the  fret-space  transit  time  from  the 
transmitting  to  the  receiving  antenna 

Convolution  of  (  5 .3  I  with  the  5i  r  I  terms  of  (  5  3  >  and  (541 
lusi  reproduce  the  function  rnirl  ((involution  w  ilh  the  other 
terms  may  become  significant  only  after  a  time  of  order  7 
Accordingly,  if  the  incident  field  has  a  time  duration  small 
compared  with  7,  i  e  ,  its  spatial  length  is  small  compared  with 
the  in-between  antennas  distance  /.  then  the  tune  dependence 
of  the  transmitted  field  is  simply  given  by  the  time  convolu¬ 
tion  of  the  incident  signal  and  the  function  70(r>  This  trans¬ 
mitted  field  is  the  same  as  would  be  obtained  for  the  case  of 
plane  wave  excitation  Accordingly,  the  result  is  obtained  that 
the  finite  distance  between  antennas  plays  no  significant  role 
if  the  incident  waveform  is  sufficiently  short  in  lime  For  in¬ 


stance,  if  the  incident  signal  is  a  pulse  of  unit  amplitude  and 
time  duration  T  ,  then 


(5  5  I 


\'r  -  T’ 


where  t*  =  /  -  ( 7/ 1  )  and  is  the  retarded  time  A  qualitative 
sketch  of  (5.5)  and  l5.bl  is  given  in  Fig.  4  for  T’  >  3r j.  When 
T’  <  3tj,  then  the  transmitted  field  is  just  given  by  (5. 3 1 
times  7 

VI  CONCH 'SIONS  AND  PRACTICAL  CONS1DF  RATIONS 

We  have  considered  the  problem  of  the  transmission  of 
steady-siate  and  transient  electromagnetic  waves  through  a 
slab  An  analytical  solution  has  been  obtained  for  the  case  of 
a  linear  homogeneous  isotropic  highly  conducting  infinite 
slab  excited  by  collinea-  electric  or  magnetic  dipoles.  The 
transmitted  :  components  of  the  field  are  expressed  as  the 
product  isteady-state  easel  or  the  convolution  (transient  easel 
of  the  corresponding  incident  field  components  and  a  two- 
term  factor  In  the  frequency  domain  the  first  term  of  this 
lactor  (see  (5  111  is  exactly  the  transmission  coefficient  of  a 
plane  wave  normally  incident  on  the  slab.  The  second  term 
takes  into  account  the  finite  distance  between  the  transmuting 
and  receiving  antennas  and  becomes  significant  only  when  this 
distance  is  ol  the  order  of.  or  sma,.er  than,  the  tree-space 
wavelength  (steady-state  easel  or  the  spatial  length  ot  ihe  inci¬ 
dent  pulse  (transient  easel  It  is  therefore  possible  to  obtain 
plane  wave  excitation  results  even  when  sources  (and  receivers) 
are  located  at  finite  distances.  For  this  all  that  is  needed  is  ihe 
proper  choice  of  distance  between  antennas 

It  is  certainly  true  that  these  results  have  been  obtained 
under  the  conditions  that  the  transmitting  antenna  is  a  dipole 
oriented  normal  to  the  slab,  that  the  transmitted  field  is  com¬ 
puted  along  the  axial  direction  of  the  dipole,  and  that  onlv  the 
z  components  of  the  field  are  used  in  the  comparison  How¬ 
ever,  we  believe  that  oar  analysis  has  a  more  general  validity 
For  instance,  results  of  the  colhnear  configuration  can  easily 
he  exlended  to  transmitted  field  points  off  the  axis  We  should 
only  suhstituie 

J  o(£\'P2  +  P2  ~~  -PP  cos  t>)  (b  I  ) 

for  J0{ fpl  in  (3  131  Then  expansion  (301  of  the  Bessel  fun  - - 
tion  <  t  !  j  and  iniegation  in  <p  gives 

f- it{p.  ! )  ~  u)Ck  J  ( 1  -  u2  )  t{u )  J  q[k  p\f~]  —  u2  ) 

*T 

•  exp  <  -  /k  lu)  du  1 1  2  i 
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which  is  the  generalization  of  (3.7)  to  the  case  p  #  0.  Then 
dF,‘(p,  1)1  dp  =  0  for  p  =  0,  which  implies  that  the  results  of  our 
analysis  are  certainly  also  valid  in  the  neighborhood  of  the 
axis.  Furthermore,  use  of  Maxwell's  equations,  with  (6.2)  as 
the  longitudinal  field,  shows  that  the  same  is  true  for  trans¬ 
verse  fields.  In  this  extension,  however,  the  simplifying  as¬ 
sumptions  used  in  the  body  of  this  paper  should  be  checked 
again.  Should  further  study  show  that  the  above  considera¬ 
tions  can  be  extended  to  more  complicated  geometries,  all 
simulation  studies  for  shielding  purposes  might  be  worth 
reconsidering. 

Some  few  practical  notes  are  now  in  order.  Reterence  is 
made  to  a  copper  slab  (o  =  5.8  X  107  S/m)  of  thickness  s  = 
I  mm,  so  that  t  —  1.52  X  10— 19  s  and  rj  =  70  ps.  Only  the 
plane  wave  transmission  coefficient  will  be  considered.  For 
incident  pulses  of  unit  amplitude  and  time  duration  T  <  r ). 
the  peak  of  the  transmitted  field  is  equal  to  2.7  X  10— 7  T'/s), 
therefore  linearly  decreasing  with  the  bandwidth  (~l/7">  of 
the  signal.  In  the  sinusoidal  excitation  case  the  attenuation  due 
to  the  mismatch  41  a  I  equals  that  due  to  the  damping  inside 
the  slab  material  exp(  -  la  lxr/\/2)  at  the  frequency  /  =  0.72 
MHz.  At  this  frequency  the  transmitted  field  is  equal  to  1 1  X 
10“ 12  times  the  incident  one.  At  higher  frequencies  the 
signal  is  decreasing  exponentially  with  the  square  root  of  the 
frequency. 

For  moderate  antenna  spacings  it  is  noted  that  the  transmit¬ 
ted  field  can  be  computed  using  the  plane  wave  transmission 
coefficient  only  when  the  attenuation  is  very  high  However, 
this  may  not  be  the  case  if  even  small  apertures  exist  in  the 
screen  Accordingly,  we  believe  it  is  worthwhile  to  extend  the 
analysis  presented  in  this  paper  to  other  canonical  problems, 
which  are  amenable  to  the  same  analytical  approach.  Among 
those,  we  list  the  problem  of  an  infinite  conductive  screen 
with  a  single  hole  and  that  of  a  conductive  screen  with  a  regu¬ 
lar  lattice  of  equal  small  apertures.  The  former  problem  can 
take  advantage  of  the  solution  of  a  plane  wave  diffraction  by 
apertures  in  conducting  screens  (211—1 23 ]  and,  eventually, 
of  symmetrization  procedures  ( 24 1 .  The  latter  could  make  use 
of  the  artifical  dielectric  theory  (251  properly  accommodated 
to  this  single  sheet  problem. 
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Abstract.  Mechanical  forces  and  torques  associated  with  electromagnetic  waves  impinging 
on  several  objects  are  computed.  The  incoming  radiation  can  be  linearly  or  circularly 
polarized,  thus  carrying  linear  and  angular  momenta.  The  objects  are  matched  dipoles  in 
several  configurations  and  a  metal  sphere.  Numerous  interesting  results  are  obtained. 

PACS:  02.41 


It  is  well  known  that  an  electromagnetic  wave  exerts 
forces  and  (possibly)  torques  on  charges  and  currents 
and  that  the  mathematical  description  of  these 
mechanical  actions  is  given  by  the  Lorentz  force 
equation,  which  connects  electromagnetic  and 
mechanical  phenomena.  If  an  electromagnetic  wave 
interacts  with,  say,  a  metal  body,  charges  and  currents 
are  induced,  and  hence  the  body  becomes  subject  to 
local  forces  and  torques.  When  a  spatial  summation  is 
made,  the  total  time-dependent  (body)  force  and 
(body)  torque  acting  on  the  object  are  obtained. 

It  is  usually  believed  that  these  forces  and  torques  are 
“small”  and  hence  negligible  from  the  viewpoint  of 
engineering  application.  However,  this  is  not  always 
so.  As  a  matter  of  fact,  a  strong  interest  in  this  area  is 
now  emerging.  For  instance,  a  book  recently  has 
appeared  [1]  which  gives  a  most  complete  review  of 
existing  theoretical  and  experimental  methods  for  de¬ 
termining  electromagnetic  signals  by  means  of  their 
associated  mechanical  effects.  Electromagnetic  pulses 
(EMP)  produced  by  nuclear  explosion  in  the  upper 
atmosphere  is  another  case  of  recent  interest.  Also 
satellite  applications  have  attracted  considerable  at¬ 
tention,  where  the  (relatively)  small  mass  of  the  object 
and  the  cumulative  nature  of  the  mechanical  actions 
(translation  or  rotation)  can  produce  long-range  mac¬ 
roscopic  effects.  Mechanical  effects  can  be  used,  for 
example,  to  steer  the  satellite  by  radiating  suitably 
polarized  electromagnetic  fields  from  the  satellite.  And 
even  key  operations,  as  deploying  antennas,  can  be 


*  On  leave  of  absence  from  University  of  Naples,  Italy 


performed,  in  principle,  by  injecting  transient  currents 
into  the  structures.  There  are  cases  in  which  the 
movement  of  a  metal  object  in  an  electromagnetic  field 
can  produce  additional  mechanical  stresses.  For  in¬ 
stance,  the  propellers  of  a  helicopter,  which  may  be 
roughly  modeled  as  rotating  dipoles,  can  be  further 
stressed  when  exposed  to  an  incident  electromagnetic 
field. 

Accordingly,  we  believe  that  it  is  useful  to  examine  a 
number  of  very  simple  cases  of  interaction  between 
metal  objects  and  electromagnetic  fields,  and  to  com¬ 
pute  the  resulting  mechanical  forces  and  torques.  As 
the  metal  object,  we  have  considered  the  simplest  one  - 
the  tuned  electric  dipole  -  in  several  configurations: 
stationary,  translational,  rotational,  receiving,  trans¬ 
mitting,  where  the  electric  dipole  is  isolated  or  coupled 
with  a  magnetic  dipole  as  a  turnstile  antenna  or  a 
chiral  antenna.  We  have  computed  the  body  force  and 
torque  acting  on  the  dipole  and  its  time  average,  for 
the  case  of  an  incident,  or  a  transmitted,  harmonic 
electromagnetic  field.  As  a  general  result,  the  forces  are 
always  proportional  to  the  Poynting  vector  (divided  by 
the  velocity  of  light)  times  an  “equivalent  area”,  which 
is  four  times  larger  than  the  usual  effective  area  - 
3AJ/4n  -  of  the  (tuned)  dipole.  This  clearly  implies  that 
the  mechanical  action  depends  on  the  power  scattered 
by  the  dipole.  Similarly,  the  torque  is  proportional  to 
the  Poynting  vector  (divided  by  the  angular  frequency 
of  the  wave)  times  the  same  equivalent  area. 

The  existence  of  the  torque,  however,  is  not  dependent 
only  on  the  symmetry  properties  of  the  scattering 
object  or  on  the  polarization  of  the  incident  wave.  We 


0340-3793 '800023 -0153  S01.80 


154 


G  Franceschetti  and  C  H  Papas 


Fig.  1  Receiving  electnc  dipole  excited  by  a  monochromatic  linearly 
polarized  electromagnetic  plane  wave 


have  examined  cases  in  which  the  incident  field  - 
circularly  polarized  -  carries  an  angular  momentum 
and  yet  exerts  no  torque  on  the  dipole.  Conversely,  we 
have  examined  cases  in  which  the  dipole  -  a  simple 
short  wire  -  radiates  as  a  scattered  field  a  linearly 
polarized  wave  (in  the  far-field)  and  has  torques  exert¬ 
ed  on  it.  The  full  analysis  of  a  rigorous  boundary 
value  problem  -  a  metal  sphere  in  an  incident  field  - 
shows  that  the  transient  force  has  not  always  the  same 
sense  as  the  incident  radiation;  in  other  words,  the 
sphere  is  alternately  pushed  and  pulled  in  the  direction 
of  the  incident  Poynting  vector. 

The  large  variety  of  effects  we  discovered  in  this  “case” 
analysis  seems  to  indicate  that  this  is  an  underde¬ 
veloped  field  of  very  interesting  research,  and  a  num¬ 
ber  of  general  theorems  could  probably  be  discovered, 
thus  bringing  to  a  better  understanding  of  a  wrongly 
neglected  field. 


1.  Mechanical  Forces  and  Torques  Exerted 

on  a  Receiving  Dipole 

by  an  Incident  Electromagnetic  Wave 


of  opposite  sign  appear  at  its  end  points.  In  (1)  the 
expression  in  the  denominator  is  the  radiation  resis¬ 
tance  of  the  (tuned)  dipole  and  J  the  free-space  intrinsic 
impedance. 

The  incident  magnetic  field,  H,  interacts  with  the 
induced  current,  thus  producing  a  force 

E2  3a2 

F  =  nldlxxH=  —  — -  cos2acos2a>rz  (3) 

£c  2n 

whose  time  average  is  given  by 
S  3A2 

<F>  =  -  A  cos2az,  A  =  —  ;  S  =  \EH  ,  (4) 

c  2n 


where  c  is  the  light  velocity. 

Accordingly,  the  force  is  given  by  the  radiation  pres¬ 
sure  S/c  times  the  “effective  area”  of  the  unloaded 
dipole,  which  is  four  times  larger  than  the  area  of  a 
matched  dipole 


GA2  3a2 
4rt  8n 


(5) 


In  addition  to  the  magnetic  force,  the  incident  electric 
field  E  interacts  with  the  induced  charges.  No  resulting 
force  is  produced  from  this  interaction;  however,  a 
torque 

T  =  Q£sin  ad/ cos  wt  sinwrz 

acts  on  the  dipole.  Its  time  average  is  given  by 

<T>  =0 .  (6) 

Accordingly,  no  time  average  torque  is  exerted  on  the 
dipole. 

The  conclusion  is  that  a  linearly  polarized  electromag¬ 
netic  wave  exerts  a  force  on  a  receiving  dipole  tuned  on 
the  frequency  of  the  wave,  where  the  force  is  pro¬ 
portional  to  the  incident  power  density  and  is  directed 
along  the  direction  of  propagation. 


/./.  Linearly  Polarized  Incident  Field 


We  examine  the  case  of  an  electric  dipole  of  length 
dl-4  A,  tuned  on  A,  and  excited  by  a  normally  incident 
plane  monochromatic  electromagnetic  wave  of 
angular  frequency  to  and  wavelength  A  (see  Fig.  1).  The 
incident  field  is  linearly  polarized;  and  the  electric 
vector  £  cos  ait  makes  the  angle  a  with  the  dipole 
direction. 

Due  to  the  incident  field,  a  current 


r_  £  cos  ad/ 

2»r{df2/3A2 

is  excited  along  the  dipole ;  and  charges 
„  / 

Q  =  —  sin  wt 
w 


U) 


(2) 


1.2.  Circularly  Polarized  Incident  Field 

Let  us  now  consider  a  circularly  polarized  plane 
electromagnetic  wave  which  can  be  conveniently  de¬ 
scribed  in  the  frequency  domain  by  the  phasors : 

E=£(x  +  iy);  H=^(y-ix)  (7) 

incident  on  the  (tuned)  dipole.  Phasors  of  currents  and 
charges  induced  on  the  dipole  are  the  same  as  those 
given  by  (1)  and  (2)  with  cos  a  *  1  and  with  the  time 
variation  suppressed. 

The  average  force  acting  on  the  dipole  is  still  that 
produced  by  the  interaction  of  the  induced  current 
with  the  y-component  of  the  magnetic  field  (the 
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x -component  is  ineffective);  its  time  average  value  is 
given  by 

<F>=i-/4i;  S  =  EH  .  (8) 

2c 

Note  that  now  the  time-average  Poynting  vector,  i.e,, 
the  incident  power  density  is  doubled  compared  to  the 
linearly  polarized  case,  and  the  average  force  is  pro¬ 
portional  to  one  half  of  the  total  incident  power. 

A  torque  is  produced  due  to  the  interaction  of  the 
y-component  of  the  electric  field  and  the  induced 
charges,  hence 

T  =  QiEAlz ,  Q=—  (9) 

ICO 

This  torque  has  the  same  sense  as  the  polarization ;  its 
time-average  value  is  given  by 

<“» 

and  therefore  is  proportional  to  a  torque  density  -  one 
half  of  the  power  density  divided  by  the  angular 
frequency  of  the  wave  -  times  the  effective  area  of  the 
dipole. 

The  conclusion  is  that  a  circularly  polarized  elec¬ 
tromagnetic  wave  exerts  forces  and  torques  on  a 
receiving  tuned  dipole,  in  the  direction  of  propagation 
and  in  the  sense  of  polarization,  respectively ;  forces  as 
well  as  torques  are  proportional  to  one  half  of  the 
incident  power  density. 


2.  Mechanical  Forces  and  Torques  Associated 
with  a  Transmitting  Dipole 


2.1.  Linearly  Polarized  Dipole 


We  examine  the  case  of  an  electric  dipole  of  moment  p, 
directed  along  the  x-axis  of  a  reference  Cartesian 
system,  as  depicted  in  Fig.  2.  The  electromagnetic  fields 
associated  with  the  dipole  are  given,  in  the  frequency 
domain,  by 


£,  = 


exp(iicr) 

4ttr 


E.  =  -  cos 
c 


E4  =  ~  sin  <)> 


(H-) 


exp(iM-) 

4nr 


„  ■  ,/■  l\exp(i*r) 

H,  =  -  nop  sin  <b  ik  —  — - - 

\  H  4w 


H,=  -iW<kos*(iK-l)e-^Kr) 


(11) 


wherein  k  —  ui/c  and  is  the  free  space  propagation 
constant  and  reference  is  made  to  the  spherical  coor¬ 
dinate  system  of  the  figure. 

The  (time  average)  force  density  transmitted  by  the 
dipole  across  the  surface  r  =  const  is  given  by  [2] 


</>  -  ie(|£.l2  +  l£J2  -  |£J2)f  +  iMIH»l2  + 1  H.\2)r 
Ipl2 


-l 


*  (4 nr)2 


2(1  -sin20cos2d>) 


sm 2  61  cos  2<p 


+  -j-y(l  —  5sin20cos2d>) 
k  r 


(12) 


Expression  (12)  is  invariant  under  the  transformation 
—  0,  4>-*n  —  <f>.  Therefore,  the  (time  average)  total 
force  transmitted  across  the  half-sphere  r  =  const,  zSO 
is  equal  and  opposite  to  the  total  force  transmitted 
across  the  other  half-sphere  r  =  const,  z  §  0  The  former 
is  given  by 


2*  nil 


<F>  =  z  J  d<{>  J  d0r2  sin  0<f>  ■  f 
0  0 


.1^1 

2  c 


1 

2  kV 


(13) 


wherein  PT  is  the  total  power  transmitted  by  the 
dipole.  Note  that  the  component  of  <F>  parallel  to  the 
plane  z  =  0  is  zero.  Accordingly,  the  (time  average) 
total  force  transmitted  across  a  half-sphere  surface  is 
directed  along  z  and  changes  sign  at  xr-1.  The 
intensity  of  this  force  is  proportional  to  the  power 
radiated  by  the  dipole  in  the  half  space  z  >0,  Pr  2, 
divided  by  the  velocity  of  light. 
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The  i time  average*  torque  density  transmitted  by  the 
dipole  across  the  surface  r  =  const  is  given  by  [3] 

:  t>  =  -  Jr  Re  {tErf  «E'  +  nH,r  x  H*J 

_  _  1  .pp  -<v  3 

2  £  (4;rr)2 

( 6  sin  9  cos  <t>  sin  <p  +  4>  cos  8  cos  <t>  sin  <j>) 


Re 


i  + 


M- 


(14) 


The  (time  average)  total  torque  <T>  transmitted  across 
the  half-sphere  is  equal  to  zero.  Accordingly,  a  linearly 
polarized  electnc  dipole  transmits  no  torque  and  only 
forces.  The  dipole  is  in  equilibrium  with  respect  to  the 
transmitted  forces. 


2.3  Circularly  Polarized  Dipole 

Let  us  now  consider  two  orthogonal  transmitting 
electric  dipoles  such  that  the  total  dipole  moment  is 
given,  in  the  frequency  domain,  by 

p  =  p(i-riy).  (15) 


With  reference  to  the  spherical  system  of  coordinates 
of  Fig  2.  the  field  components  associated  with  the 
dipole  are  given  by 


E, 


-sinfl(cos^  +  isin 

£ 


2i*c 


2  \  exp(iicr) 
r1)  dnr 


E. 


-  cos  0  (cos  <t>  +  i  sin  <j>) 


exp(ixr) 

4nr 


r  P,  •  i*  I  2\  exp(iw) 

Es  =  -(-sin d>  +  icosd>)  — - j  +  K2  - ; - 

*  £  \  r  rl  )  4  nr 

Hg  =  -itup(sin^-icos^))(iK— 


(16) 


Ht  =  —  icupcos0(cos04-isin  <p)  iK 


exp(  — ixr) 
4  nr 


The  (time-average)  force  density  transmitted  by  the 
dipole  across  the  surface  r= const  is  now  given  by 


t!e!1 

W  4  (1SF  £ 


2(1  +  cos20) - j-Tsin20 


+  -t-jU  +  cos20-4sin20) 
k  r 


r  . 


(IT) 


Also  in  this  case,  and  due  to  the  symmetry  of  (17),  the 
total  (time  average)  force  transmitted  across  the  half 
sphere  r*  const,  zgO,  is  balanced  by  an  equal  and 
opposite  force  transmitted  across  the  complementary 
half  sphere  r  =  const,  z^O;  and  hence  the  dipole  is  at 
rest.  The  total  (time  average)  force  transmitted  for  z SO 


is  given  by 


<F>  = 


2c  l*  *c2r2 


2 


z  . 


118) 


Again,  the  total  (time  average)  force  transmitted  is 
along  the  z-direction  and  is  proportional  to  one  half  of 
the  total  transmitted  power  divided  by  the  velocity  of 
light.  Note  that  for  the  antenna  under  consideration 
the  radiation  resistance  is  doubled  compared  to  that  of 
the  single  dipole. 

The  (time  average)  torque  density  transmitted  by  the 
dipole  is  given  by 


(19) 


The  total  torque  transmitted  across  the  half-sphere 
r  =  const,  z SO  is  readily  computed  as 

<T>  =  t .  (20) 

2oi 


The  torque  transmitted  across  the  complementary 
half-sphere  is  the  same.  Accordingly,  the  dipole  is  not 
in  equilibrium  as  far  as  torques  are  concerned. 


3.  Mechanical  Forces  and  Torques  Associated 
with  a  Moving  Dipole 

3.1.  Translating  Dipole 

We  examine  the  case  of  an  electric  (tuned)  dipole, 
excited  by  an  electromagnetic  plane  wave,  as  depicted 
in  Fig.  1.  We  are  taking  the  angle  a  to  be  zero  for 
convenience,  and  we  assume  the  dipole  to  move  with  a 
constant  velocity  v  =  -  ui  along  the  negative  z  axis. 

It  is  evident  that  the  force  acting  on  the  dipole  is  given 
by  the  same  expression  (4)  with  a =0,  provided  that  the 
Poynting  vector  S  is  computed  in  a  reference  frame 
fixed  to  the  dipole. 

If  E\  H'  are  the  electromagnetic  fields  in  this  frame,  we 
have 

E'  =(E  +  v  x  pH) 

H'  =  (H-VX£E)^’ 

wherein  p  *  v/c.  Then,  after  a  few  computations,  we  get 
the  following  expression  for  the  Poynting  vector  in  the 
primed  frame  (fixed  to  the  dipole) 


l 


(21) 


(22) 
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This,  however,  does  nor  imply  that  the  average  force 
exerted  upon  the  dipole  is  increased  by  the  scalar 
factor  appearing  in  (22)  in  comparison  to  the  case  of 
the  dipole  at  rest  As  a  matter  of  fact,  the  incident  wave 
in  the  primed  reference  frame  has  a  Doppler  shifted 
frequency 

l2J> 

Accordingly,  the  effective  area  of  the  dipole  -  see  (4)  - 
takes  the  new  value 

A-Arr,  ,24> 

When  the  dipole  is  retuned  on  the  new  frequency,  the 
force  acting  upon  it  remains  constant. 


Fig  3.  Receiving  rotating  electnc  dipole  excited  by  a  linearly 
polarized  electromagnetic  field 


3.2.  Rotating  Dipole 

We  examine  now  the  case  of  an  electric  (tuned)  dipole, 
excited  by  an  electromagnetic  plane  wave,  as  depicted 
in  Fig.  3  The  dipole  is  rotating  with  constant  angular 
velocity  (2  along  the  z-axis.  so  that  the  dipole  makes 
with  respect  to  the  x-axis  an  angle  <t>  given  by 

tp  =  dt .  (25) 


We  can  compute  the  forces  and  torques  acting  on  the 
routing  dipole  in  the  following  manner.  We  transform 
the  pertinent  results  of  Sect.  1  from  the  inertial  labo¬ 
ratory  frame  of  Sect.  1  to  the  non-inertial  rest  frame  of 
the  routing  dipole  and  then  neglect  manifestly  small 
correction  terms.  Usually  the  procedure  of  transform¬ 
ing  from  inertial  frames  to  non-inertial  frames  or  vice 
versa  is  a  long  calculation,  but  in  the  case  of  present 
interest  simplifying  approximations  can  be  introduced 
and  the  procedure  can  accordingly  be  shortened.  In 
particular,  we  are  interested  in  the  components 


Ee  =  Ecost/> ;  E^=  -Esintfi ; 


(26) 


which,  upon  transformation  to  the  rotating  reference 
frame,  are  given  by  [4] 


£.  =  £.;  e;= 


h.  ■  »•  = 
V^f2'  * 


1/1-0 


2  ’ 


wherein  the  parameter 
0  =  fl<?/c 


(27) 


(28) 


ukes  into  account  the  relativistic  correction.  However, 
when  the  fields  are  evaluated  at  the  rotating  elemen- 
ury  dipole,  we  have  d(-» 0,  e-»0  and  0-*O. 
Accordingly,  the  force  acting  on  the  dipole  is  given  by 


F  = 


E2 

—  A  cos2<(>i  . 
sc 


(29) 


In  order  to  consider  the  time  average  of  (29),  it  should 
be  remarked  that  the  fields  are  changing  in  time  with 
an  angular  frequency  to  and  so  the  quantity  to  be 
averaged  is  cos2  tot  cos2  dt.  If  the  average  is  made  with 
respect  to  the  period  2n/to  of  the  incident  wave,  then, 
after  some  lengthy  algebra,  we  get 


<F>  **  2c  A‘ 


sin 


2nd 


1  + 


to 


2nd 

to 


cos  2(2t  + 


2nd 

to 


d/co 

l-(d/co)2 


sin  (2tuf)sin 


2dt  + 


2nd , 
lo  I 


(d/to)2 
1  —  (d/to)2 


cos  (2a>t)  cos 


2dt  - 


2nd 


(30) 


If  d/to4  1,  the  term  in  the  curly  brackets  is  approxi¬ 
mately  equal  to 


1  +  cos(2flr) ; 


(31) 


and,  after  another  time  average  with  respect  to  the 
period  2n/d,  we  recover  the  result  (8).  When  d/to  =  1. 
the  term  in  the  curly  brackets  equals  3/2. 

To  evaluate  the  torque  acting  on  the  dipole,  we  first 
compute  the  charges  (of  opposite  sign)  at  the  endpoints 
of  the  dipole,  hence 

C-2^V3Fi“S“COs0''" 

EAl  to  sin  wt  cos  dt -d  cos  tot  sin  dt 
"  2n(Al2/3i2  to2-d2  '  (32) 

Then,  the  (instanuneous)  torque  is  given  by 
T=Q£,d/f 

E2  .  _  to  sin  wt  cos  dt -d  cos  wt  sin  dt 

=  —  Azcostof  sinflt - 5 — — i - . 

{  to2-  d2 
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Fig.  4  la)  and  lb)  Two  chiral  dipoles 


Fig  5  Geometry  of  interaction  between  the  magnetic  Field  and  the 
loop  of  the  chiral  dipole 


By  similarly  averaging  over  the  period  2 n/<o,  we  get 


<T>=—  A 


w/(2 


2w  1  -  ((2/air 


l+((2/<o)J 


InQ 


sin 


l- 


U) 


2nd 

w 


cos 


2nQ\ 

(2°'+~ ) 


,  tt  cos  2cot  cos  2(2t  v - . 

1  -{Q/wr  \  w  ) 


212/ai 

L-((2/oi)J 


sin  2ai(  sin 


(34) 


When  (2/oi<l,  the  term  in  the  curly  brackets  is 
approximately  equal  to 


1  -cos(2(2t)  (35) 

and,  after  another  time  average,  we  recover  the  result 
( 10)  weighted  by  the  factor  (2/oi  4 1.  When  (2/oi  - 1,  the 
term  in  the  curly  brackets  is  proportional  to 
[1  -((2/oi)]J;  therefore,  the  average  torque  is  zero. 


4.  Mechanical  Forces  and  Torques  Exerted 

on  a  Receiving  Chiral  Dipole 

by  an  Incident  Electromagnetic  Wave 

4.1.  Linearly  Polarized  Incident  Field 


Let  us  consider  two  chiral  dipoles  [5]  as  depicted  in 
Fig.  4.  We  will  define  as  right-handed  the  chiral  dipole 
“ft”,  while  the  chiral  dipole  “a"  is  defined  as  left-handed. 
In  addition,  the  dipoles  will  be  defined  as  balanced  if 
the  following  relation  exists  between  the  length  Al  of 
the  (electric)  dipole  and  the  are  na2  of  the  (magnetic) 
dipole 


Al  =  xNna2 


2  n2Na2 
A 


(36) 


wherein  N  is  the  number  of  turns  of  the  loop  Under 
these  circumstances,  the  far  field  radiated  by  the  chiral 
dipoles  is  everywhere  circularly  polarized,  and  the 
radiation  resistance  is  twice  as  large  as  that  of  the 
single  electric  dipole. 

We  examine  now  the  case  of  a  linearly  polarized 
electromagnetic  held  described,  in  frequency  domain, 
by 

E  =  Ex,  H=fy  |37) 


and  incident  upon  the  (tuned)  dipoles  Each  dipole  will 
be  subject  to  a  z-directed  force  which  can  be  computed 
as  in  Sect.  1,  with  the  only  difference  being  the  new 
radiation  resistance.  Hence 


<F>=~/4z,  (38) 

2C 

wherein  A  is  given  by  (4),  as  usual.  Equation  (38) 
coincides  with  (8). 

Now,  at  variance  with  the  single  dipole,  the  current 
circulating  in  the  loop. 


EAl 

4/rCdiJ/3Ai 


(39) 


will  interact  with  the  incident  magnetic  field,  thus 
producing  an  (instantaneous)  torque  (see  Fig.  S) 


T  = 


E2  M  Na 2  ^/pz 

Ta  di 


J  ±  4>  sin  4>d<t> 
o 


(40) 


Here  use  has  been  made  of  (36)  and  the  upper  and 
lower  signs  refer  to  the  right-handed  and  left-handed 
dipoles  respectively.  Accordingly,  the  (time  average) 
torque  is  given  by 

<T>*±ils/l  (4l) 

which  is  identical  to  (10). 
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4.2 .  Circularly  Polarized  Incident  Field 


We  examine  now  the  case  of  an  incident  wave  cir¬ 
cularly  polarized  along  the  z -direction,  i.e.,  in  complex 
notation, 

E  =  £(i  +  iy);  H-= E(y-ix)  .  (42) 

v 

This  wave  is  incident  upon  a  balanced  (tuned)  chiral 
dipole,  either  right-handed  or  left-handed.  The  open 
circuit  voltage  induced  by  the  electric  field  (along  the 
electric  dipole)  is  given  by  EAl.  Similarly,  the  open 
circuit  voltage  induced  by  the  magnetic  field  in  the 
loop  is  given  by 

IN  ^(  — io»)-£d/ ,  (43) 

upon  use  of  the  balance  condition  (36).  Then,  exam¬ 
ination  of  Fig.  4  clearly  shows  that  the  two  voltages 
add  for  the  right-handed  chiral  dipole,  while  they 
cancel  out  for  the  left-handed  one.  Accordingly,  one  of 
the  two  dipoles  is  “invisible”  to  the  incident  circularly 
polarized  wave,  and  no  force  or  torque  is  acting  upon 
it.  On  the  contrary,  the  current 


EAl 

2<d/J/322 


(44) 


is  induced  on  the  other  dipole.  This  current  will 
interact  with  the  incident  magnetic  field  thus  produc¬ 
ing  an  average  force 

<F>  =  Az  (45) 


S.  A  Rigorous  Boundary  Value  Solution. 
Mechanical  Action  on  a  Metal  Sphere  Produced 
by  an  Electromagnetic  Pulse 


We  examine  the  case  of  a  metal  sphere  excited  by  a 
plane  electromagnetic  wave  as  depicted  in  Fig  6  For  a 
steady-state  incident  field,  the  total  field  components 
on  the  spherical  surface  r=a  are  given  by  [6] 

E  cos  ip  “  Pi  (cos  0) 

E'  (ica)2  £"(  1)12  '[Kah^dca)]' 


//0=  |  T,(e)sin<f> ;  W*=  t  T2(0)cos<p  , 


(49) 


equal  to  (8).  Note  that  the  forces  due  to  the  interaction 
of  the  magnetic  field  with  the  current  along  the  loop 
have  a  zero  spatial  average. 

Now,  we  compute  the  torque  acting  on  the  dipole.  This 
is  produced  in  part  by  the  y-component  of  the  mag¬ 
netic  field  interacting  with  the  current  along  the  loop. 
Hence 

Ti=z^—A  (46) 

and,  in  part,  by  the  y -component  of  the  electric  field 
interacting  with  the  charges  at  the  end-points  of  the 
wire.  Hence 

E2 

Tj=  z  —  A  .  (47) 

C  oj 

Accordingly,  the  average  total  torque  is  given  by 
<T>-i^A  (48) 


wherein 


T(9)=  —  T  iir"2n+  if  1  p-(cosg) 

’  KaV"  Mn+lUtKahiVa)]'  sin0 


i  <fP‘(cos0)l 
_____  —  | 


r.w-i.iMr-'iill 


1 


[Kah^a)]' 
Pi  (cos  0)1 


(50) 
dPi(cosO) 

— a _ 


d6 


[icahj/^Ka)]'  sin  0 


and  a  prime  implies  the  derivative  operation  with 
respect  to  the  argument1. 

Let  £( to)  be  the  spectral  distribution  of  the  incident 
field,  i.e.. 


£»E0/(co)exp(i»cz) ,  (51) 

wherein  £0  is  a  convenient  normalization  constant, 
and  t  «a/c.  Let  us  further  introduce  the  following 


which  is  twice  as  large  as  (10). 


1  Note  ft  misprint  enor  in  th<  expression  for  £,  in  [6]. 
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time-dependent  functions 


,  *,*  ^(tu)exp(-ia>t)(  — i)"'1' 

"  J.  (a)T)[ojTh,.ll(<Dt)] 

*5  ^(a))exp(-ia)t)(-  l)** 
{wT)[unh["(wty]i 

7.(0=  ;  1  f.U)dt  . 

I  -  « 


1 

-da) 


2 

—  da) 


(52) 


Then,  the  transient  field  at  the  sphere  surface  r  =  a  is 
given  by 


£,  =  -£0cosd>£«(2n+  l)7.(MP,!(cos®) 
1 


u  _E0 sin<t>*  2n  +  l  P.1  (cos 6) 

"•  -  ;  V  I"" 


=  £p  cos  0  *  2h-f1 

*  C  V'nin+i) 


7.(0 


dPj  (cos  6) 
d® 


+  9.(0 


P'  (cos  ®) 
sin® 


(53) 


Transient  surface  density  charges 


acting  on  the  sphere  has  only  a  z-component  given  by 

F.  =  i  7  Z.  L»  (2n  +  D(2m  + 1)/.(07.(0 

C  1  1 


* 

•  Jcos®sin®P,|(cos®)P1l1(cos®)d®,  (59) 

o 

wherein  S  =  E\/ (. 

By  repeated  use  of  recurrence  and  orthogonality  re¬ 
lations  [7],  it  can  be  shown  that  the  integral  appearing 
in  (59)  is  different  from  zero  only  for  m  =  n  +  1  and 
m  =  n  -  1.  Then,  rearranging  the  terms  of  the  series,  we 
get 

F,  =  z|4naJ|^,n(n+l)(n  +  2)7jt)7..i(0  •  (60) 

c  i 

For  the  second  density  force  we  proceed  similarly.  The 
x  and  y  components  of  the  (space-average)  total  force 
acting  upon  the  sphere  are  zero,  and  only  the 
z-component  is  different  from  zero.  This  last  com¬ 
ponent  is  given  by 


r  .  S  2  sr-  ^  2n+ 1  2m  + 1 

L.L.- 


,  n(n  +  1)  m(m+  1) 


(/./»  +  9,9 J 


?  /  Pi  Pi  dP\  dP'm\ 
(sin  0  sin  0  +  dO  dO  ) 


cos® sin  Odd 


*  /  p>  dp' 

(fd)m  +  fmSj  j  + 


sin  6  d6  I 


e4  =  eEr 


(54)  •  cos® sin  Odd 


(61) 


and  surface  density  currents 

J,  =  —H4  ;  J.  =  H,  (55) 

are  excited  on  the  sphere  surface.  The  electric  field 
will  interact  with  the  charges,  thus  producing  a  force 
density 

f,  =  eE,2f  (56) 

and  the  magnetic  field  will  interact  with  the  currents, 
thus  producing  a  force  density 

f + 

For  the  first  force  density  we  have 

f,  =  E£oCOS24>f£,  Yj*  (2n+  l)(2m+  l)J,U)Jjt) 

i  i 

Pi  (COS®)  Pi(C0S®) 
sin  ®  sin  0 

This  force  density  has  a  cos3<£  dependence  for  its 
x -component;  and  a  cos2 #  sin  4>  dependence  for  its 
y-component.  Therefore,  the  space-average  of  these 
two  components  is  zero.  This  does  not  happen  for  the 
z-component  of  (58),  whose  ^-dependence  is  of  type 
cos2$.  Accordingly,  the  (space  average)  total  force 


(57) 


(58) 


The  integral  which  appears  in  (61)  can  be  evaluated 
upon  repeated  use  of  recurrence  relationship  between 
associated  Legendre  polynomials  and  of  their  defining 
differential  equation  [8],  In  particular,  the  second 
integral  is  different  from  zero  only  for  n  =  m;  and  the 
first  one  only  for  m  =  n+  1  and  m  =  n—  1.  Then,  upon 
rearranging  the  terms  of  the  series,  we  get 


c  1 


n(n  +  2) 


n+  1 


(/./.♦  1+9.9  n+  l) 


2n+  1 
+  n(n+  1) 


(62) 


Accordingly,  the  (space-average)  total  force  acting  on 
the  sphere  is  given  by 


P  ® 

F=  —  z~dna2Y- 
c  , 


nfn+2) 

n+1 


(/./.+ 1 +9.9  .+ 1) 


+ 


2n+  1 
n(n+  1) 


fjS, ~  "("+!)("  +  2)7,7,. , 


(63) 


An  interesting  case  is  that  of  a  bandlimited  signal  such 
that  (oma  <  1 ,  wm  being  the  maximum  angular  fre¬ 
quency.  Then,  using  the  series  expansion 


xhj,1  *(x) 


exp(ix)  (2n) ! 
r*  1  n ! 


t 

(2ilr 


(64) 


1 
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it  follows  that,  to  the  lowest  order  in  zd/dt 


£o9.(f)-(-)"'l(2rr 


"! 2  y-»E 


(2n)!  dt"-' 


(65) 


where,  now,  £  =  £(r)  and  is  the  time-varying  incident 
field.  Accordingly,  at  the  lowest  order  again 


F  =  —  i  tt-  Ana2  2r£(t*)  ~ 

CC  1*1“.*  “  0 

4  na2  <fS  I 


(66) 


where  S(r)  is  the  instantaneous  value  of  the  incident 
Poynting  vector  at  z  =  0.  Note  that,  at  the  lowest 
considered  order,  the  time  average  of  F  is  equal  to  zero 
for  a  sinusoidal  incident  field,  since  the  sphere  does  not 
scatter  any  real  power.  For  a  Gaussian  pulse 


£(:.  t)  =  £0  exp 


U  +  z/c)2 
2T 2 


(67) 


FU>  =  2 


2zt* 
T 2 


(68) 


It  is  seen  that  the  sphere  is  first  pushed  (r*  <0)  in  the 
sense  of  propagation  of  the  incident  wave;  and  then 
pulled  (f*  >0). 


6.  Conclusions 

We  have  considered  the  mechanical  effects  that  an 
electromagnetic  wave  can  have  on  various  types  of 
targets.  Using  nothing  more  than  the  Lorentz  force, 
given  by  F  =  pE  + J  x  jiH.  we  have  found  that 

a)  a  linearly  polarized  wave  falling  on  a  tuned  electric¬ 
receiving  dipole  will  exert  a  body  foce  -  see  (4)  -  on  the 
dipole  in  the  direction  of  incidence, 

b)  a  circularly  polarized  wave  falling  on  such  a  dipole 
will  exert  not  only  a  body  force  but  also  a  torque  -  see 
(10). 

We  also  found  that 

c)  for  a  fixed  transmitting  dipole  surrounded  by  an 
absorbing  spherical  surface  of  radius  r,  there  is  no 
recoil  force  on  the  dipole  but  there  is  a  force  on  the 
absorbing  spherical  surface  -  see  (13)  -  which  tends  to 
make  it  implode  when  icr<  I  and  explode  when  kt>  1. 
Moreover, 

d)  when  the  fixed  transmitting  antenna  happens  to  be 
a  turnstile,  there  is  a  recoil  torque  on  the  antenna,  and 
a  torque  on  the  absorbing  spherical  surface  -  see  (26), 


and  again  there  are  the  implosive  and  explosive  forces 
as  before. 

in  addition,  we  have  studied  the  problem  of  a  receiving 
dipole  that  is  either  moving  at  constant  velocity  or 
simply  rotating.  We  have  found  that 

e)  when  the  tuned  dipole  is  in  pure  translational 
motion,  the  force  exerted  on  the  dipole  by  the  incident 
signal  is  equal  to  that  on  the  tuned  dipole  at  rest. 

f)  When  the  dipole  is  only  rotating,  the  incident  signal 
exerts  also  a  torque  on  the  spinning  dipole  -  see  (34) 
For  a  balanced  chiral  dipole  receiving  antenna  in  the 
field  of  a  Linearly  polarized  incident  wave,  we  deduced 
that 

g)  the  body  force  exerted  on  the  antenna  is  inde¬ 
pendent  of  the  handedness  of  the  antenna  -  see  (38), 
whereas  the  torque  changes  sign  with  a  change  of 
handedness.  If  such  an  antenna  is  exposed  to  a  cir¬ 
cularly  polarized  incident  wave,  then 

h)  the  antenna  is  completely  free  of  body  force  and 
torque  for  one  handedness,  but  not  for  the  other  -  see 
(45)  and  (48).  That  is,  with  the  proper  handedness  the 
antenna  is  "invisible"  at  the  “balanced”  frequency. 
Finally,  the  push-pull  effect  has  been  examined  for  the 
case  of  a  pulse  falling  on  a  spherical  target  It  has  been 
found  that 

i)  the  leading  term  of  body  force  excited  on  the  sphere 
by  the  incident  pulse  is  proportional  to  dS  cl,  provid¬ 
ed  that  the  effective  bandwidth  of  the  pulse  multiplied 
by  the  characteristic  time  of  the  target  is  less  than 
unity. 
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